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SUMMARY 


This report describes some detailed trip distribution analyses 
of the 1971 census journey to work data for the fifteen Ontario census 
areas other than the Toronto CMA. Earlier studies of the commuting 
patterns in these fifteen Ontario urban areas had demonstrated that the 
commonly used gravity model is capable of explaining only the broad pat- 
terns of interaction between home and work and that significant over - 
and under - estimation residuals existed between the observed and simu- 
lated trip interchanges. A number of potential opportunities for im- 
proving the gravity model were identified in the earlier study and the 
principal objective of the current study was to explore the impacts of 
these modifications to the gravity model structure. 


The cell entries in observed home to work linkages matrices 
reflect two types of effect and these are the rate of interaction between 
a pair of zones and the absolute magnitudes of the labour force and the 
employment in those zones. The essential hypothesis of the gravity 
model is that the rate of interaction between zones is a function of the 
spatial separation of zones. In order to understand better the spatial 
interaction patterns in the Ontario census areas it is essential to 
extract from home to work linkages matrices the pure spatial interaction 
effects that are independett of zone size effects. 


Bi-proportional »alancing techniques have been used to trans- 
form the home to work linkages matrices so that the entries in each row 
and in each column sum to a constant magnitude. The effect of this 
operation is to produce matrices in which the entries represent the in- 
teraction magnitudes that would occur between zones of equal size, or 
the pure interaction effects between home and work zones. This balancing 
process is straight forward in that the observed matrix row entries are 
each scaled by the same proportion so that the row totals are all equal 
to an arbitrary constant. The columns of this adjusted matrix are then 
summed and scaled so that the column totals are all equal to the arbi- 
trary constant. This row by column adjustment process is continued 
iteratively until convergence is obtained with each row total and each 
column total being equal to the arbitrary constant. If the arbitrary 
constant were equal to one then the cell entries in the adjusted matrix 
could be interpreted as the unit probabilities of interaction between 
zones. 


A clustering technique was then used to group zones on the 
basis of similarities in their destination vectors of the bi-proportionally 
adjusted matrix. The technique used is Ward's method which uses a 
measure of the Euclidean distance between destination vectors to cluster 
census tracts. The procedure begins with the n census tracts in a 
census area and groups them into a sequence of clusters until all census 
tracts are merged into a single cluster. The heterogeneity of the census 
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tracts within a cluster may be represented by the within cluster sum 
of squared differences in the destination vectors of the census tracts. 
The sum of the sum of squared differences over all clusters provides 
an indication of the overall heterogeneity in the commuting patterns 
from the residential zones! at each stage of the clustering process. 
The error sum of squares increased slowly during the first stages of 
the clustering process and began to increase sharply when residential 
zones with quite different destination vectors are forced into the 
same cluster. The hierarchical structure of the clustered residential 
zones in a particular census area is conveniently represented graphi- 
cally by a dendrogram. The dendrograms show the sequence in which 
zones cluster along with the error sums of squares magnitudes at which 
they cluster. 


Dendrograms are presented for the fifteen Ontario census 
areas and commuting sub-regions are identified in each area on the 
basis of the cluster analysis. A detailed review of these commuting 
sub-regions suggests that there are five broad categories of deter- 
minants of these regions and these are: (i) multi-community compo- 
sition of a census area; (ii) topographic and man-made features; 
(iii) the time sequence of development of an area; (iv) socio-economic 
factors; and (v) the domination of large employment concentrations. 


Thunder Bay, Kitchener, St. Catharines and Hamilton census 
areas provide examples of the impact of municipal boundaries on com- 
mutersheds. In almost all cases of clustering at the four or five 
cluster level the cluster boundaries are coincident with municipal 
boundaries. The best examples of the impacts on commuting patterns 
of topographic features are provided by London, Hamilton and Ottawa. 
Timing of development influences are noted in almost all of the census 
areas particularly those that grew rapidly during the 1966-1971 period. 
The strongest determinant of commuting patterns appears to be the 
socio-economic characteristics of the households and in many cases the 
other determinants mentioned above co-vary with the socio-economic 
influences. More formal analyses of the determinants of commuting 
patterns could be undertaken by performing cluster analyses of house- 
hold characteristics and exploring the extent to which these cluster 
boundaries are coincident with the commutershed boundaries. 


Previous studies of the gravity model have shown that the 
commonly used goodness-of-fit statistics are inadequate for assessing 
the extent to which gravity model calculated trip matrices agree with 
the observed trip matrices. The effectiveness of a range of statistics 
is explored by a series of simulation experiments for eight of the 
larger Ontario census areas. Simulated trip matrices were calculated 
from the observed trip matrices by multiplying each of the observed 
cell entries by a randomly generated error with a specified range. A 
sequence of simulated trip matrices were prepared for six error ranges 
from ten percent to one hundred and fifty percent and magnitudes of 
the goodness-of-fit statistics were calculated for each of the six 
error ranges. 
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Detailed analyses of the changes in the magnitudes of the 
goodness-of-fit statistics with increasing induced error magnitude have 
shown that the so-called phi statistic of information theory is the most 
appropriate statistic to use. This statistic increased linearly with 
increasing error magnitude and resulted in a consistent ordering of 
the eight census areas across all census areas. While the phi statistic 
magnitude increases with increasing numbers of trips at a constant error 
magnitude a normalized phi magnitude was found to be roughly constant 
across all census area sizes. Accordingly the phi statistic divided 
by the number of trip interchanges in an area has been adopted as the 
most appropriate goodness of fit measure for evaluating changes to the 
gravity mddel structure. 


A large number of variations in the gravity model have been 
tested for the Ontario census areas. A very flexible gravity model 
computer program has been developed which allows a wide range of gravity 
model structures to be tested and very exhaustive goodness-of-fit 
statistics to be calculated. In most cases gravity models have been 
calibrated for the eight largest Ontario census areas ranging in size 
from Thunder Bay to Ottawa. 


Production constrained, attraction constrained and doubly 
constrained versions of the gravity model have been calibrated using 
road network distances in kilometres as the surrogate for the genera- 
lized travel costs between zones. This basic set of models was cali- 
brated in order to provide the datum against which changes in the 
gravity model might be assessed. The calibration criterion used for 
all model types is the minimization of the sum of the absolute dif- 
ferences in the ordinates of the observed and simulated trip length 
frequency distributions. In most cases this criterion also satisfied 
the more normally used criterion of equality of observed and simulated 
area-wide mean trip lengths. 


Using the normalized phi statistic as the goodness of fit 
measure of the calibrated models suggested that these models produced 
trip tables that were of the same statistical quality as those pro- 
duced by the random error generator with an error range of 75 to 100 
percent. A detailed comparison of the model types showed that the 
doubly constrained model was consistently superior to the production 
and attraction constrained versions of the gravity model with the 
attraction constrained version having the worst goodness of fit. These 
analyses also showed that the goodness of fit deteriorated significantly 
with increasing urban area size. While the doubly constrained model 
provided the best goodness of fit there are extreme difficulties in 
providing any reasonable interpretation of the travel deterrence para- 
meter and balancing factor magnitudes. 


Doubly constrained models using network travel times rather 
than network distances have been calibrated for the Kitchener, Hamilton 
and Ottawa census areas and compared with those developed in terms of 
network distances. The phi-statistics showed that these models em- 
ploying travel times were marginally inferior to those calibrated in 
terms of network distances. Analyses of the zone-specific goodness of 
fit statistics showed that in the Ottawa CMA the time-based models 
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improved the goodness-of-fit of some of the outlying residential 
census tracts but that this was not a general trend across all census 
areas. 


Inter-zonal generalized travel costs were also created from 
various combinations of network distances and network times for the 
Kitchener CMA. Analyses of the goodness of fit statistics showed that 
doubly constrained models calibrated using these generalized cost func- 
tions are superior to those using travel times but marginally inferior 
to the models calibrated using network distances. 


Detailed comparisons of the trip interchange residuals for 
the Ottawa CMA and a partitioning of the phi-statistic into over - and 
under - estimation residuals by intra-zonal and inter-zonal trip inter- 
changes showed that the major source of error is the under-estimation 
of trip interchange magnitudes and that the superior behaviour of the 
doubly constrained model is due to a reduction in the under-estimation 
residuals. However this is not a completely unbiased interpretation 
of the residuals since earlier simulation studies had shown that the 
phi-statistic is more sensitive to under-estimation errors than to 
over-estimation errors. 


A number of stratified gravity models have been calibrated 
using all three basic forms of the gravity model. The three groups of 
stratified gravity models calibrated are: (i) multi-parameter gravity 
models with separate travel deterrence function parameters estimated 
for each census tract; (ii) sub-region specific models in which 
separate travel deterrence‘ parameter magnitudes are calibrated for 
from four to eight commuting sub-regions depending on the specific 
‘area; and (iii) a set of models calibrated for separate socio-economic 
groups. 


The multi-parameter production constrained models all ex- 
hibited superior goodness of fit characteristics to the single parameter 
production constrained models but were not as good as the single para- 
meter doubly constrained models. The largest improvements occurred in 
the multi-community census areas of Kitchener and St. Catharines. The 
multi-parameter attraction constrained models were inferior to the single 
parameter versions except for Thunder Bay, Kitchener and St. Catharines. 
The multi-parameter doubly constrained models all exhibited superior 
goodness of fit characteristics to the single parameter versions. 


The cluster analyses of the census areas formed the basis 
for the identification of the calibration sub-regions where these 
ranged from three in Thunder Bay and Kitchener to eight in the Ottawa 
census area. The sub-region specific production constrained models 
are superior to the single parameter versions but are inferior to both 
the single and multi-parameter production constrained models. An 
analysis of the sub-region specific doubly constrained models showed 
that they performed marginally better than the single parameter version 
but are inferior to the multi-parameter models. 


The most useful models stratified by socio-economic char- 
acteristics are those developed for non - car owners and car - owners. 
Detailed studies for the Kitchener area showed that doubly constrained 
versions of these stratified models perform at a superior level to the 
Single parameter doubly constrained models. 


Trip generation rates have been established for each of the 
census areas for households stratified by tenancy status and by period 
of residence. Owner occupied dwelling units have significantly higher 
numbers in the labour force than the rented dwelling units. The number 
of home to work linkages from households with periods of residency 
greater than five years is larger than for households with shorter 
periods of residency. 


Detailed regression analyses have been conducted in order to 
establish transit captive and non-captive rates for all census areas. 
Statistically significant prediction equations have been developed in 
terms of the dwelling unit composition of census tracts and there are 
significant differences in the rates between census areas. Transit 
captive generation rates from single attached dwelling units are from 
two to three times greater than from single detached dwelling units. 
Some of the variation in generation rates is due to the variation in 
the labour force participation rates between census areas and a second 
set of regression equations has been developed from the data factored 
to a constant labour force participation rate across all census areas. 
Pooled regression equations were also developed for the entire set of 
census areas. The pooled equations show that for average labour force 
participation rates the captive generation rate is 2.3 times larger for 
attached dwelling units than for detached dwelling units and that the 
nor-captive rate is about 3 times larger for detached dwelling units 
than for attached dwelling units. Significant relationships between 
captive and non-captive trip attraction and employment sector type 
could not be established. 
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INTRODUCTION 


The practical importance of the journey to work data collected 
in the 1971 Census of Canada has been recognized for many years by the 
Ministry of Transportation and Communications of Ontario. This organ- 
ization has been responsible for the use of census information in a 
number of transport model building and research projects conducted both 
within the Ministry and by external agencies. This report describes 
the results of a research project sponsored at the University of Waterloo 
by the Ministry of Transportation and Communications which has been con- 
cerned with a more detailed analysis of the spatial interaction patterns 
in Ontario census areas observed in the 1971 census. This project re- 
presents an extension to an earlier research project conducted during 
1976-1977 which was also concerned with the analysis of the census data. 

The primary objective of these two research projects was to 
explore the use of the census journey to work data in the development 
of generalized transport models that might be used to synthesize future 
work trip demands without the collection of comprehensive special 
purpose surveys. A second objective of these studies was to develop a 
better understanding of spatial interaction patterns across a wide range 
of urban areas and the factors that influence these interaction patterns. 

The research project completed in January, 1977 showed that 
regression equations with very high statistical qualities could be 
developed from the census data for estimating the amount of labour 
force in each census tract from a knowledge of the dwelling unit com- 


position of census tracts. Various forms of the gravity model were 


calibrated using the census journey to work data. While these cali- 

brated models explained much of the observed spatial interaction 

between census tracts large positive and negative trip interchange 
residuals were observed in all urban areas. It was suggested that 

a variety of factors influenced these residuals and these included 

the timing of development, socio-economic factors and the nature of 

the deterrence function embodied in the gravity model. 
In an attempt to develop better trip distribution models 

a number of specific objectives were set for the research project 

described in this report and these included: 

1. To use bi-proportional matrix balancing techniques to calculate 
normalized trip linkage magnitudes between each pair of origin- 
destination zones in order to separate interaction effects from 
zone size effects 

2. To use clustering techniques to identify groups of census tracts 
with similar interaction behaviour in order that factors influencing 
spatial interaction may be isolated and incorporated in the gravity 
model structures 

3. To develop a more flexible gravity model computer program that 
allows alternative model types incorporating different deterrence 
functions to be calibrated and model errors to be examined in 


detail 


4. To examine various statistical criteria that may be used to assess 
the goodness of fit of trip distribution models 


5. To estimate the parameter magnitudes of alternative forms of the 
gravity model which have been calibrated by geographic sub-region, 
socio-economic group and which incorporate various travel deterrence 
functions and travel deterrence measures. 


CHAPTER 1 


THE CENSUS JOURNEY TO WORK DATA 


The characteristics of the 1971 journey to work data and 
the urban road networks which have been coded to supplement the 
spatial interaction data are described in detail in ref. [1]. How- 
ever, it is worthwhile recalling the general characteristics of the 
data base used in the research project described in this report. 

The basic spatial unit to which the 1971 census journey to 
work data are coded is the census tract where the average residential 
census tract contains about 5,000 people. In 1971 the number of census 
tracts varied from 14 in Guelph (population = 62,660) to 452 in the 
Toronto CMA (population = 2,628,000). ‘The coded journey to work data 
represent approximately an 11 percent sample of the journey to work by 
households. The contiguous urbanized portions of a census area are 
covered by the official census tract system and the analyses described 
in this report deal with the interactions between official census tracts: 
Home to work linkages also occur between official census tracts and the 
Surrounding census districts and with other census areas but these inter- 
actions are excluded from the analyses. In most of the Ontario census 
areas the home to work linkages between official census tracts within a 
census area account for about 80 percent of the total interaction. 

Road networks have been coded at roughly the arterial/sub- 
arterial level in terms of road distances. In addition, link travel 


speed data have been added for selected urban areas. 


1.1 The Ontario Census Areas 

At the 1971 census there were 16 census areas in Ontario which 
ranged in population size from Guelph to Toronto. This report is con- 
cerned with the journey to work in 15 of these census areas with the 
Toronto CMA being excluded from the analysis. Table 1 shows the pop- 
ulations of the 15 census areas in 1971 along with the number of census 
tracts in each area. This table shows that six of the census areas had 
populations of less than 100,000, three had populations between 100,000 
and 200,000, four had populations between 200,000 and 300,000 and two 
had populations of 500,000 or greater. 

In addition to the variations in population size the census 
areas also varied widely in economic base, spatial character and geo- 
graphic setting. In 1971 three broad industry groups accounted for about 
70 percent of the jobs in the Ontario census areas and these were manu- 
facturing, trade and services. The principal differences between the 
economic bases of the census areas were in the primary (mining), manu- 
facturing, transportation and public administration sectors. In Sudbury, 
26.5 percent of the jobs were in mining while Brantford, Oshawa and 
Kitchener had strong manufacturing industry bases. Employment in the 
transportation asim was well above the average in Thunder Bay while 
16.3 percent of the labour force in Ottawa was in public administration. 

These variations in the employment bases of the Ontario census 
areas were reflected in the spatial distributions of employment. The 
geographic settings of the census areas also have a strong influence on 
their spatial characters. Eight of the Ontario census areas have heavy 


concentrations of employment in two or three census tracts while in four 


TABLE 1. 1971 Census of Canada Census Metropolitan Areas and 
Census Agglomerations for Ontario 

C(M)A 1971 Population Number of 
Number Total C(M)A Oriicial’ tis CTs CDs 
Po: Guelph 62,660 60,087 14 Lg 
2 Peterborough 63,530 Byers Bilt 16 10 
3 Sarnia 78,445 60,417 18 6 
4 Brantford 80,285 70,904 16 16 
5 Sault Ste. Marie S270 805352 igs) Z 
6 Kingston alate Y hs 59 S047 18 8 
7 Thunder Bay £22,095 Ihe} sya SLI 25 2 
8 Oshawa 12075320 im ea a 22 10 
9 Sudbury 1S 5425 94,624 28 5 
10 Kitchener 226, 850 212,819 45 ai 
11 Windsor 258,645 211,494 56 5 
tz London 286,010 PAROS ie RMP Sy 5 
LS St. Catharines 303,430 Pele 2 Oe 53 Lt 
14 Hamilton 493,925 437,554 109 Wg 
its) Ottawa 602,510 506,316 120 19 


of the census areas residential areas are located away from the job 
centres creating longer average trip lengths. In several census areas 
these effects combine to produce long average trip lengths that are 


unrelated simply to population size. 


1.2 Previous Trip Generation Analyses 


Ref. [1] describes in detail the regression analyses conducted 
in the 1976-1977 research project which were directed towards the develop- 
ment of equations for predicting the labour force in each census tract 
given a knowledge of the amount of residential activity in a tract. It 
must be remembered that in this report the term trip generation really 
means the number in the labour force in each zone rather than the number 
of trips. The best prediction equations were developed in terms of the 
dwelling unit composition of census tracts and these prediction equations 
are summarized in Table 2. This table shows that the equations are of 
very high statistical quality and that the partial regression coefficients 
are quite stable across urban areas. The consistency and explanatory 
power of these equations would suggest that work trip production studies 
be standardized in terms of dwelling unit composition. Trip attraction 
analyses could not be performed with the aggregate data since independent 


measures of employment by census tract were not available. 


1.3 Previous trip Distribution Analyses 


Three forms of the tradional gravity model were calibrated 
from the census data and these were a production constrained model incor- 
porating an inverse power deterrence function and a production constrained 


model with attraction trip end balancing incorporating inverse power and 


TABLE a2; Census Tract Labour Force Prediction 
Equations Using Dwelling Unit Composition 


Census Area a by b, R 
Guelph oi) te Ifa Ps Oe 
Peterborough 0 1.544 0.909 0.98 
Sarnia =63 12526 0.983 OSE: 
Brantford 26 15.506 0.925 G29 
Sault Ste. Marie -51 Moye!) LoD ae, 
Kingston -28 1.669 IDR ahfe tlh Oro2 
Thunder Bay 45 1.438 Os919 0.96 
Oshawa -26 1.524 fe bo 0.96 
Sudbury 24 152) 1.394 (ORS 
Kitchener =87 ee Ow 1.310 0523 
Windsor -—]7 eos OF 919 O97 
London 197 1.474 L026 0.97 
St. Catharines 98 L435 OU79'6 Ogg 
Hamilton 13 Hells 1.083 Be, 
Ottawa anh 1.660 te 212 0.96 


CT Labour Force = a + by: CT Single Detached Dwelling Units 


a5 b,- CT Single Attached Dwelling Units 
and Apartments 


negative exponential deterrence functions. While there were some differ- 
ences in the overall performance of the various model types all of the 
calibrated models produced large residuals between simulated and observed 
trip interchange magnitudes. 

Figure 1 illustrates the spatial distributions of under- and 
over-estimation residuals for the production constrained version of the 
gravity model in London. The trip interchange residuals plotted in 
Figure 1 are equal to the estimated trip interchange magnitudes minus 
the observed trip interchange magnitudes. Inspection of the residuals 
presented in Figure 1 By pe Vath most of the under-estimation residuals 
are associated with non-CBD employment locations while almost all of the 
over-estimation residuals are associated with the CBD. Employment in 
London in the central census tract is very much higher than in the 
other census tracts and in the production constrained model trips tend 
to be diverted to the central zone from competing employment locations. 

Figure 2 illustrates the spatial distributions of under- and 
over-estimation trip interchange residuals for the doubly constrained 
form of the gravity model calibrated for London. This diagram demon- 
strates that the strong over-estimation of trips to the central area is 
reduced significantly by the use of the attraction trip end balancing 
procedure. The balancing procedure reduced the attractivity of the central 
area and increased the attractivity of the non-CBD employment locations. 
Figure 2 illustrates that while trip end balancing reduced the majority 
of over-estimation residuals it had little impact on the under-estimation 
residuals. Similar effects to those illustrated in Figures 1 and 2 for 


London were observed for other Ontario census areas. The size of the 
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FIGURE 1. Production Constrained Gravity Model Trip 
Interchange Residuals for London 


- 10 - 


LONDON 


KILOMETRES 


OVEREST IMAT IONS 


LONDON 
KILOMETRES 


400 


LOWER MASK= 75 
UPPER MASK= 10000 
UNDEREST IMAT IGNS 


FIGURE 2. Doubly Constrained Gravity Model Trip Inter- 
change Residuals for London 
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residuals indicated that significant modifications to the gravity model 
structure are required if spatial interaction magnitudes are to be esti- 


mated more realistically. 


1.4 Additional Analysis Opportunities 


In the earlier research project the gravity models were cali- 
brated using the aggregate home to work matrix. This approach assumes 
implicitly that the interaction patterns observed in 1971 are influenced 
only by the spatial distributions of production and attraction trip ends 
and the spatial separation of zones. It has been suggested previously 
that this simple hypothesis is inadequate and that other factors have 
an important influence on home to work linkages. 

The gravity models traditionally calibrated have a Single deter- 
rence function which implies that the response to travel time of residents 
in all parts of an urban area is similar. Figure 3 shows the home to work 
linkages greater than 250 between census tracts in the Kitchener census 
area. This diagram illustrates that commuting tends to be relatively self- 
contained within each of the four local municipalities and that the character- 
istic trip lengths tend to vary between areas. This suggests that one 
potential improvement to the gravity model might be to estimate behavioural 
parameters for the separate sub-regions of the area. 

A second characteristic of urban areas that is not reflected in 
the traditional form of the gravity model are the differences in commuting 
behaviour between various socio-economic groups. Clearly the interaction 
spaces of people without access to a car for the journey to work will be 
restricted. Figure 4 shows the spatial distributions of productions and 


attractions for car-owners and non-car owners in the Kitchener census area. 
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This diagram illustrates that households containing non-car owning members 
of the labour force and jobs attracting nor-car owners are not distributed 
uniformly throughout the area but tend to be ee around the central 
area of Kitchener. A second potential improvement to the gravity model 
structure is to stratify it by socio-economic group in order to reflect 
potential differences in commuting behaviour between groups. 

A third factor influencing spatial interaction patterns is the 
staging of development of an urban area. The location decisions of individ- 
uals at a particular time will reflect the spatial distributions of house- 
hold and employment opportunities that exist at that particular time. It 
would be expected that households which have located in recent times would 
tend to be concentrated on the periphery of the urbanized area where new 
residential sub-divisions are developing. In addition one would expect 
that these new residential areas would be linked strongly with those areas 
of new employment growth. These types of effects are illustrated in 
Figures 5 and 6 for the Kitchener census area. Figure 5 shows the home to 
work linkages for households which have been located for six years or more, 
that is, prior to 1966. In the Kitchener area most of these linkages are 
between households located in the middle and inner suburbs and the inner 
industrial areas with some linkages to jobs in the industrial park in the 
southern sector of Kitchener. Figure 6 illustrates that the spatial link- 
ages of households which have been located for five years or less have a 
different character with the bulk of the trips originating from the outer 
residential areas. Many of these trips are to the University of Waterloo 
and the Kitchener CBD both of which grew rapidly during the 1966-1971 period. 


A third opportunity for improving the predictive powers of the gravity model 
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is to recognize that home to work linkages develop in a sequential way 
and not instantaneously as implied by the traditional gravity model. 
Clearly, each of the three types of effects mentioned previously are not 
mutually exclusive and there are some overlaps between behaviour, socio-— 


economic characteristics and the timing of development. 
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CUAPTERS2 


EMPIRICAL ANALYSES OF SPATIAL INTERACTION 


The previous chapter has identified a number of opportunities 
for improving the ability of the gravity model to explain observed 
spatial interaction between home and work. The gravity model estimates 
spatial interaction in terms of the product of the production and 
attraction trip end magnitudes and the travel deterrence between zones. 
Some of the difficulties of improving the gravity model are related to 
this confounding of the zone size effects and the interaction effects. 

The cell entries in observed trip matrices reflect these two 
types of effect, the rate of interaction between a pair of zones and 
the labour force and employment magnitudes of the zones. In order to 
understand better the spatial interaction patterns it is essential to 
extract from the home to work linkages matrices the pure spatial inter- 
action effects that are independent of zone size effects. It is these 
rates of spatial interaction that trip distribution models attempt to 
explain in terms .of the travel deterrence function. 

This chapter describes some empirical studies of the journey 
to work matrices obtained in the 1971 census for the fifteen Ontario 
census areas. The analysis strategy consisted of two basic steps which 
are: (i) bi-proportional balancing of the trip matrices; and (ii) the 
clustering of zones in terms of similarities in destination characteristics. 

The bi-proportional balancing techniques transform the home to 


work linkages matrices so that the entries in each row and in each colum 


make yere 


Sum to a constant magnitude. The effect of this operation is to produce 
matrices in which the entries represent the interaction magnitudes that 
Would occur between zones of equal size, or the pure interaction effects 
between home and work zones. The clustering technique compares the des- 
tination vectors for each origin zone and groups zones together that have 
the most similar destination vectors. The clustering technique groups 


all zones into a single cluster using a step by step approach. 


area: Bi-Proportional Matrix Balancing 


It has been suggested above that bi-proportional balancing 
techniques may be used to remove the effects of zone size from trip 
matrices leaving only the information on the rate of spatial inter- 
action between zones. If this is not done then zones with large amounts 
of employment will dominate any analyses and may mask important inter- 
action effects between smaller zones. Bi-proportional matrix balancing 
techniques have been developed for use in a number of fields with much 
of the original work being conducted in connection with input-output 
tables [2]. The Furness method [3] of origin-destination table develop- 
ment is a closely related procedure. 

The entries in bi-proportionally adjusted trip matrices re- 
present the interactions that would occur if all the zonal productions 
and attractions in an urban area were equal. The row/column constant 
used in matrix balancing can have any arbitrary Magnitude. For example, 
if it were equal to 1.0 then each cell entry could be interpreted as 
the probability of interaction between a known origin zone and a part— 


icular destination zone when all destination attractivities have the 


same magnitude. Normalized trip matrices therefore reflect only spatial 
separation effects rather than both spatial separation and structural 
effects. 

The calculation of bi-proportional matrices is straight forward. 
Trip matrix row entries are each scaled so that the row totals are all 
equal to the arbitrary constant. The columns of this new matrix are 
then summed and scaled so that the colum totals are all equal to the 
arbitrary constant. This i ge is continued iteratively until con- 
vergence is obtained with all row and column totals being equal to the 
arbitrary constant. 

It is useful to examine the properties of the transformed 
matrices in relation to the original trip matrix. The bi-proportional 
balancing techniques preserves the interaction rate structure of the 
original matrix but the cell entries in any row or colum are not all 
adjusted by a factor with a constant magnitude. However it has been 
shown [2] that the products of the final row and column adjustment 
factor magnitudes have unique magnitudes for each cell. This is, the 
interaction magnitudes calculated by the bi-proportional balancing 
techniques are unique. 

It is also useful to examine the changes in the trip length 
structure that might be expected from the transformation of the original 
trip matrix. The mean trip length of any origin or destination zone is 
calculated simply from the weighted sum of the trip distances to each 
of the other zones in the system. It has been noted previously that 
the bi-proportionalized matrix retains the interaction rate structure 


but not the interaction magnitude structure of the original matrix. 


As a result the mean trip lengths for any origin or destination zone for 
the bi-proportionalized matrix may be expected to change since the row 
or column weight composition used will have altered. Rows which do not 
have strong concentrations of destinations at one or two zones may be 
expected to exhibit little change in the zone mean trip length. On the 
other hand zones with strong concentrations of destinations in a few 
zones might be expected to exhibit Significant shifts in the mean trp 
length. The weights between zones with small trip end magnitudes will 
increase and a greater influence on the zone mean trip lengths of the 
bi-proportionalized matrices. The bi-proportional adjustment process 
will produce more uniform zone mean trip lengths than Originally observed. 

Table 3 summarizes the area-wide trip lengths for both the 
original trip matrix and the balanced matrix for the fifteen Ontario 
census areas. This table shows that in most cases the area-wide mean 
trip lengths of the balanced matrices are shorter than the observed. 

The largest shifts occur in those census areas with unusual trip linkage 
characteristics such as Sarnia, Kingston, Sudbury, Hamilton and Ottawa. 
In these areas a few employment zones tend to dominate the commuting 
Structure, 

Figures 7 to 9 provide more detailed comparisons of the ad- 
justed and observed trip matrices in which the adjusted mean trip length 
for each zone is plotted against the observed mean trip length. Figures 
7 and 8 provide comparisons for the origin zones in six census areas while 
Figure 9 shows comparisons for the destination zones of Hamilton and 


Ottawa. 
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TABLE 3. Comparison of Area-Wide Observed and Balanced Trip Matrix 
Mean Trip Lengths 


Observed Balanced 
Mean Mean 
Trip Length Trip Length 
km km 
Guelph | 3.76 Bf7 
Peterborough 36 3 3.36 
Sarnia 6,33 5.44 
Brantford 4.20 2B ail 
Sault Ste. Marie 4.59 4.55 
Kingston 5.68 4.32 
Thunder Bay ene 5.65 
Oshawa 4.10 4.10 
Sudbury 8.63 109m 
Kitchener 4.80 5.03 
Windsor 7,90 FARE! 
London 6.57 ele Ae 
St. Catharines 6.50 Sis28) 
Hamilton 8.05 7.10 


Ottawa 7.63 6.88 
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FIGURE 7. Adjusted vs Observed Origin Zone Mean Trip 
Lengths. for Kitchener and London 
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FIGURE 8. Adjusted vs Observed Origin Zone Mean Trip 
Lengths for Hamilton and Ottawa 
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Figures 7 to 8 illustrate the tendency for shorter trip lengths 
to become longer and the longer trip lengths shorter. The largest shifts 
in the mean trip lengths are associated with zones located on the peri- 
phery or the smaller-sized zones. Figure 9 illustrates that the mean 
trip lengths for the destination zones tend to be reduced by the balancing 
procedure. It should be remembered that the aim of the bi-proportional 
balancing procedure is not to preserve the trip length structure but the 


interaction structure of a community. 


2.2 Clustering Techniques 


The clustering technique used in this investigation is a pro- 
cedure known as Ward's method which is a hierarchical agglomerative 
technique [4]. This technique begins with the n census tracts in a 
census area and groups them into a series of clusters with aggregation 
continuing until all of the census tracts are grouped into one cluster. 
With this technique the analyst is not constrained to some arbitrary 
pre-determined number of clusters. 

The clustering technique uses a measure of the Euclidean 
distance between the destination-zone similarities of two census tracts 


which is defined in the following way: 
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where d.. = the Euclidean distance between the destination 

= characteristics of the two census tracts i and j 

k = the destination zone under consideration for a 


particular origin zone (i or j) 


ee 


* 
Toy = the normalized number of trips between origin zone i 
and destination zone k 
n = the total number of census tracts 


hates, ac is a measure of the Similarity in the destination vectors 
of the two particular origin zones being considered for a potential 
merger. If the destination vectors were identical then ar =a) ave 

The lower the magnitude of a the greater the Similarity between the 
two census tracts. 

The clustering technique proceeds by first calculating an upper 
triangular matrix of the similarities in the destination vectors of each 
pair of census tracts. The pair of census tracts with the minimum ve 
magnitude are then merged into a cluster. The entries in the matrix of 
Similarity are then re-calculated to reflect this merging of two census 
tracts and the procedure is continued until all census tracts are merged 
into one cluster. The logic of the clustering procedure is illustrated 
in Figure 10. 

The calculations of the new entries in the similarity matrix 
do not involve a complete re-calculation of the entries in the similarity 
matrix since not all entries are affected. Suppose two census tracts 
i and j join to form a new group ij , usually labelled with the 
smallest census tract number. The similarity between this new cluster 


and another cluster, or census tract, is given by: 
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FIGURE 10. Steps in the Clustering Procedure 
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n = the number of census tracts in the cluster identified by 
the sub-script 


Equation (2) indicates that the new measure of similarity is calculated 
from the weighted sum of the Euclidean distances between the pairs of 
clusters used in the previous step. The larger n, is relative to ae 
then the greater di will be weighted in the calculation of the new 


similarity measure d - Further, the larger the value of n then 


h?ij 
the smaller the weighting of oh which is the distance between the two 
census tracts being merged. 

The distance value calculated by equation (1) may be interpreted 
as the within cluster sum of squared differences between the destination 
vectors of two census tracts. As more census tracts are included in a 
cluster then the distance value as calculated by equation (2) represents 
the within cluster sum of Squared differences between a cluster centroid 
and the location of the individual census tracts in the cluster. The 
clustering procedure groups census tracts so as to minimize the increase 
in the sum of squared differences ecross the entire set of census tracts. 
One would expect that the first Stages in the clustering hierarchy would 


consist of a number of pairwise census tract clusters each with a small 


increase in the sum of squared differences. Large increases in the aggre- 
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gate sum of squared differences might be expected as clusters are grouped 


to form eventually a single cluster. 


2.2.1 An Example 


Figure 11(a) shows a simple five zone home to work matrix of 
the type available from the census. This matrix may be transformed into 
the doubly-balanced matrix shown in Figure 11(b) by alternatively scaling 
the row and column entries so that they sum to 20, an arbitrary constant. 

The entries in the matrix of Figure 11(b) may be used along 
with equation (1) to produce the upper triangular matrix of similarity 
measures shown in the upper part of Figure 12. For example: 
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and so on. 
Inspection of the initial similarity matrix of Figure 12 shows 
that the minimum similarity coefficient is d,. with a magnitude of 2.49. 
The first step in the clustering procedure, then, is to group zones 3 
and 5 to form a new zone 3' , where the marginal increase in the error 
be updated to reflect this initial clustering of zones 3 and 5 and this 
UB? aDAle aw! 


be used to calculate these new similarities: 


requires the calculation of d . Equation (2) may 
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FIGURE 11. Doubly Balanced Home to Work Linkages Matrix 
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The revised Similarity matrix which reflects the merging of zones 3 and 


5 is shown in Figure 12. It should be re-called that q represents 


Lu 
the Similarity between zone 1 and the new zone 3' formed from the clustering 
of zones 3 and 5. 

Inspection of the revised Similarity matrix shows the minimum 
distance is ds, =)2,.05 and. the second step in the clustering procedure 
is to merge zones ] and 2 to form a new zone 1' with a marginal increase 
in error sum of Squares of 2.65. A new Similarity matrix has to be pro- 
duced from the calculation of REL and div, and this matrix is 
shown in the lower part of Figure 12. The minimum distance measure in 
this revised matrix is doy, which means that the third clustering 
step involves merging 3' and 4 to form a new zone 3" with a marginal 
increase in the error sum of squares of 4.5, Once again a new similarity 


matrix has to be calculated which has one CHLELY SORE d =97938% “The 


Lis 


final fusion into one cluster occurs with a marginal increase in the 


error sum of squares of 7.38. The clustering sequence for this example 
is zones 3 and 5, then zones 1 and 2, then merged 3 and 5 with 4 and 
finally merged 1 and 2 with 3, 4 and 5. 

Figure 13 shows the dendrogram for this particular example 
where this dendogram has been plotted using the convention of the com- 
puter program used in this study. The clusters are arranged as far as 
possible in terms of increasing zone number magnitude. The horizontal 
lines representing the fusion of two census tracts, or clusters, is 
drawn at the magnitude of the marginal increase in the sum of squares 
at which this union occurs where the magnitude is shown as the ordinate. 
That is, zones 3 and 5 are clustered first at a magnitude of 2.49; zones 1 
and 2 are then clustered at a magnitude of 2.65, and so on. The lower 
part of Figure 13 traces the marginal increase in the error sum of squares 
as the number of clusters decreases to 1. This diagram is simply a re- 
statement of the vertical spacings of the horizontal lines shown in the 


upper part of the diagram. 


2.3 Dendrograms for Ontario Census Areas 


The appendix to this chapter shows the dendrogram constructed 
for each of the fifteen Ontario census areas. These dendrograms have 
been plotted using the convention introduced in Figure 13. The dendro- 
grams trace the clustering hierarchy and identify the magnitudes of the 
error sums of squares at which merging occurs. The dendrograms for the 
Guelph and Brantford census areas are illustrated in Figure 14 while the 
increases in the error sums of squares with decreasing numbers of clusters 


are shown in Figure 15. 


He ee 


DENDROGRAM OF CLUSTER 


HIERARCHY 

feel 
Oc 
<q 
=) 
‘ej 
Y) 
LL 
oO 
= 
=D) 
ae) 
. 2 3 5 4 
Oo 
x ZONE NUMBER 
LJ 
a 
= 8 
LJ 
(op) 
<—{ 
LJ 
jon 
© 
= 
=| 
Ss 
aA 
ac 
<q 
=> 

O 5 4 3 2 l 


NUMBER OF CLUSTERS 


FIGURE 13. Dendrogram of Cluster Hierarchy and Associated 
Marginal Increase in Error Sum of Squares 


- 36 - 


Orr VN ye Feng 
gi 


Oo 
—-O-— OF WU Om Ow = 


BRANTF GRO 


FIGURE 14. Dendrograms for Guelph and Brantford 


- 37 - 


y GUELPH 


ite) te) 5 O 
NUMBER OF CLUSTERS 


IN ERROR SUM OF SQUARES 
Oo 


recalls 

of BRANTFORD 
LJ 

Gc 

O 

Zz 

=) 

<q 

= 

© 

Oo 

P< § 

= 

15 10 5 O 


NUMBER OF CLUSTERS 


FIGURE 15. Marginal Increase in Error Sum of Squares 
with Decreasing Numbers of Census Tracts 
for Guelph and Brantford 


Eee sake 


Inspection of the dendrogram for Guelph in Figure 14 shows that 
the first pair of census tracts to cluster are 4 and 13, followed by 7 
and 8, 5 and 12, and 9 and 11. As the clustering process continues 10 
clusters with the already clustered 9 and 11. This dendrogram shows that 
clustering continues without major jumps in the error sum of squares. 
With three clusters the census tract composition of each cluster is: 

Tsk, 1007-34... lan 5 andalZ 
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The clustering process in the Brantford region has slightly 
different characteristics than for Guelph. While the initial stages of 
the clustering are similar the error sum of squares increases sharply 
during the final stages of the clustering process. The dendrogram for 
Brantford suggests the following census tract composition for three 
clusters: 
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It is interesting to note that census tract 11 in Brantford is almost 
exclusively an ehatetat zone and that it clusters with the remaining 
census tracts only at the last stage. It should also be noted from 
Figures 14 and 15 that the increase in the error sum of squares is 
large when clusters I and II are merged and that a similarly large 
increase occurs when census tract 110 is merged with the other members 
of cluster II in the previous stage. Census tract 110 embraces the 


the Paris area which is part of the Brantford census area and its 


forced merger with the census tracts in the northwestern part of Brantford 
produces the relatively large increase in the error sum of Squares at the 
anti-penultimate stage of clustering. 

The characteristics of the Guelph and Brantford dendrograms are 
sufficient to illustrate the general properties of the dendrograms contained 
in the appendix of this chapter. Some of the effects highlighted above are 
more extreme in the larger census areas than in the two smaller census 
areas. For example, the Kitchener census area consists of the three major 
municipalities of Waterloo, Kitchener and Cambridge which embrace relatively 
self-contained commuter sheds. When these areas are forced together large 
increases in the error sums of Squares occur. The ordinates of all of the 
dendrograms presented in the appendix are plotted to the same scale and 
differences in the clustering properties of each census area are detected 
easily. The residential zone clustering properties of each of the fifteen 
Ontario are discussed in detail in the following paragraphs. It should be 
recalled that the residential zones have been clustered on the basis of 
similarities in destination vectors and that no census tract adjacency con- 
straints have been used. 

In reviewing the commuting structure of each census area five 
broad categories of effects have been used to assist in this analysis and 
these categories are: 

Multi-Community Composition: several censSuS areas embrace a number of 
separate municipalities where there is some degree of commuting self- 
containment. 

Topographic and Man-Made Influences: many census areas are bisected by 
rivers and other topographic features as well as certain man-made features 
such as rail lines which may influence commuting patterns. 

Timing of Development: strong linkages may be observed between residential 


and employment zones that develop at about the same time even though these 
zones may be separated by some distance. 


eee se 


Socio-Economic Factors: certain dwelling unit types and employment sectors 
May attract higher than average proportions of particular socio-economic 
groups and their specific locations will condition spatial linkage patterns. 


Specific Employment Concentrations: several census areas have very strong 
concentrations of employment of a specific type in one or two census tracts 
and specific commuting patterns may be associated with these areas. 


Guelph 


Figure 16 shows the census tract grouping in Guelph for three clusters. 
Cluster I consists of the six census tracts along the eastern side of the 
area, cluster II consists of the four census tracts on the western side 

and cluster III embraces the CBD (CT 6) and the two census tracts immediately 
to the west of the CBD. The peripheral clusters I and II are linked with 
job opportunities in CT 1 (University of Guelph) and CT 9 (suburban industrial 
area) while the census tracts in cluster III are linked primarily to jobs 

in the CBD. The census tracts marked with an * in Figure 16 experienced 
significant population growth during 1966-1971 and developed strong link- 
ages with employment in CTs 1 and 9. The primary factor influencing 

spatial linkages in Guelph is the timing of development with the inner 
suburbs oriented to CBD employment and the newer outer suburbs linked 

with suburban employment. Figure 17 shows the major home to work linkages 
observed in Guelph in 1971 which also illustrates these two broad com- 
muting pattern differences. 


Peterborough 


Figure 18 shows the census tract grouping in Peterborough for three clusters. 
CT 14had a very small population in 1971 and it does not cluster naturally 
with the other Peterborough census tracts. The dendrogram in the appendix 
shows that it clusters with CT 110 which includes the Lakefield Community 
north of Peterborough. Peterborough had two major concentrations of employ- 
ment in 1971 and these were the CBD in CT 7 and the manufacturing employment 
in CT 4. Cluster I includes the older peripheral residential zones adjacent 
to the manufacturing plants in CT 4 along with CTs 11 and 12 which experienced 
sigificant population growth during 1966-1971. The principal determinants 

of this cluster are strong linkages with jobs in CT 4. In contrast the CTs 
in Cluster II have strong linkages with jobs in CT 7. The CTs in this 
cluster include the older established residential areas in CTs 5, 6, 7, 8 

and 10 as well as the newer residential areas in CTs 9 and 13. The principal 
determinant of spatial linkages in Peterborough seems to be socio-economic 
with some influences of the timing of development. Figure 19 shows the 

major home to work linkages in Peterborough. 


Saineid 


Figure 20 shows the three major clusters of census tracts in Sarnia. Cluster I 
consists only of CT 1 which contains the concentration of petro-chemical jobs 

in Sarnia and only a small number of households. Cluster II contains CTs 2,3,5, 
6, 7 (CBD) and 100 (railway jobs). These CTs tend to have strong linkages 

with jobs in CTs 7 and 100. Cluster III embraces the remaining CTs along 

the eastern boundary of Sarnia which cluster on the basis of their strong 
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FIGURE 16. Census Tract Clusters for Guelph 
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FIGURE 18. Census Tract Clusters for Peterborough 
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FIGURE 19. Major Home to Work Linkages for Peterborough 
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FIGURE 20. Census Tract Clusters for Sarnia 
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linkages with jobs in CT 1. Some separation of the CTs in cluster III may 
be detected from the Sarnia dendrogram in the appendix where differences are 
due primarily to linkages with jobs in CTs 7 and 100. The principal deter- 
minants of spatial linkages in Sarnia are socio-economic factors along with 
the strong concentration of jobs in the petro-chemical industry. Figure 21 
illustrates the major trip linkages in Sarnia and these highlight the 
importance of the three major employment concentrations. 


Brantford 


Figure 22 illustrates the census tract composition of three clusters for 
Brantford where cluster III contains only CT 11 a major manufacturing zone 
with few households. Figure 22 shows that Brantford divides into two 
fairly distinct commuting clusters. While the significant amounts of 
employment in the central area CTs 5 and 6 are common to both commuter 
sheds the differentiation between the commuter sheds is created by the 
linkages to jobs in CTs 1 and 2 for cluster I and linkages to jobs in 
CTs 3 and 11 for the CTs in cluster II. The dendrogram in the appendix 
illustrates that there is a significant increase in the error sum of 
squares when clusters I and II join. The dendrogram shows that this 
large increase in error sum of squares is due to CT 110 (Paris) joining 
cluster II just prior to its union with cluster I. The principal deter- 
minants of the clusters in Brantford are its multi-community composition 
and the presence of the Grand River which provides part of the boundary 
between clusters I and II. Figure 23 illustrates the principal home to 
work linkages in Brantford. 


Sault Ste. Marie 


Figure 24 shows two major clusters of census tracts for Sault Ste. Marie. 
Employment in 1971 was concentrated in CT 8 (steel mill complex) and CT 6 
(CBD). Because of the magnitude of the jobs concentrated in CT 8 all re- 
sidential zones have relatively strong linkages with this CT. The CTs in 
cluster I are in the eastern part of the area and they have strong linkages 
with the CBD employment as well as the steel mill complex. CTs in cluster II 
have very strong linkages to the steel mill complex. The dendogram in the 
appendix shows that there is more heterogeneity in the destination character- 
istics of CTs in cluster II than in cluster I and that there is a sharp in- 
crease in the error sum of squares when clusters I and II are merged. The 
principal determinant of commuting patterns in Sault Ste. Marie is the large 
concentration of jobs in the steel mill complex with socio-economic factors 
influencing some of the linkage patterms. Figure 25 illustrates the major 
trip linkages in Sault Ste. Marie and shows the dominance of the steel mill 
oriented trip linkages. 


Kingston 


Figure 26 shows three major clusters of census tracts for the Kingston area. 
CT 13 is essentially a manufacturing employment zone which only clusters with 
the other CTs in the final stage creating a sharp increase in the error sum 

of squares. Cluster I contains the CTs in the older inner suburbs of Kingston 
which have strong linkages with Queen's University (CT 2), the hospitals (CTs 
2 and 1) and the CBD (CT 1). Clusters II and III embrace CTs that are located 
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FIGURE 22. Census Tract Clusters for Brantford 
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FIGURE 26. Census Tract Clusters for Kingston 
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to the north and south of Princess Street. While the CTs in these two 
clusters focus on different employment centres there are strong socio- 
economic differences between the two clusters of census tracts. Socio- 
economic factors are the most important determinant of home to work link- 
ages in Kingston. Topographic constraints also influence commuting link- 
ages from the residential areas in Kingston Township but these are not 
shown. Figure 27 shows the major home to work linkages in Kingston and 
illustrates the dominance of employment at Queen's University/Kingston 
General Hospital and the CBD. 


Thunder Bay 


Figure 28 shows the two major clusters of census tracts for Thunder Bay 
where these clusters embrace CTs within the former Cities of Fort William 
and Port Arthur, respectively. CTs 20 and 23 and outlying e's. 110 andel.20 
merge with the Fort William CTs in cluster I only after a significant in- 
crease in the error sum of squares. CT 20 is essentially an employment 
area and the residences in this CT exhibit an unusually long mean trip 
length. The dendrogram illustrates that destination patterns within each 
cluster are very homogeneous and that a very large increase in the error 
sum of squares occurs when clusters I and II are joined. Clearly, the 
major determinant of commuting linkages in the Thunder Bay area in 1971 
were the former municipal boundaries. Within the former Fort William area 
employment in the CBD (CT 6) and CT 3 are the principal determinants of 
commuting patterns. The self-contained nature of the commuting patterns 
is further illustrated in Figure 29 which illustrates the home to work 
linkage patterns. 


Oshawa 


Figure 30 shows the three major census tract clusters for Oshawa where the 
boundaries of clusters I and II are essentially the municipal boundaries 

of Oshawa and Whitby, respectively. The dendrogram shows that cluster III 
(Gls 105s, 110) merges with cluster II and that there is a sharp increase 
in the error sum of squares when clusters I and II join. Commuting patterns 
within Oshawa are dominated by the manufacturing employment concentrations. 
Municipal boundaries are the principal determinant of commuting patterns 
within the Oshawa census area and socio-economic factors tend to determine 
the spatial linkages patterns within Oshawa. In addition, the growth of 
housing opportunities in CT 105 north of Highway 401 has influenced the 
Spatial linkages. Figure 31 illustrates the major home to work linkages 

in Oshawa. 


Sudbury 


Figure 32 shows the Sudbury census tracts grouped into two major clusters. 
The CTs in cluster I located in the southwest sector of the area have strong 
commuting linkages with the CBD. The CTs in cluster II are located in the 
northeast and tend to have strong linkages with the mining jobs. The dendro- 
gram in the appendix identifies a third cluster of residential areas which 
are located in the townships surrounding Sudbury. A significant increase in 
the error sum of squares occurs when clusters I and II are clustered with 
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FIGURE 32. Census Tract Clusters for Sudbury 


- 60 - 


the CTs in the surrounding townships. The principal determinants of com- 
muting in Sudbury are socio-economic differences along with the locations 


of mining employment. Figure 33 shows the major trip desires in Sudbury 
and demonstrates that significant linkage magnitudes exist only to the 
employment in the CBD. 


Kitchener 


Figure 34 shows four major census clusters for the Kitchener census area 
where the cluster boundaries follow closely the former municipal boundaries 
of Waterloo, Kitchener, Preston and Galt. The dendrogram inthe appendix 
illustrates that the CTs in Preston and Galt cluster at about the same 
error sum of squares magnitude as the Kitchener CTs. A sharp increase 

in the error sum of squares occurs when the Kitchener and Waterloo census 
tracts are merged and a very high increase occurs when the Kitchener- 
Waterloo and Cambridge (Preston, Galt) clusters are merged. While manu- 
facturing employment tended to dominate the economic base of the census 
area in 1971 there were major concentrations of non-manufacturing employ- 
ment in the Kitchener CBD and at the University of Waterloo. As a result 
socio-economic factors have an important influence on commuting patterns 
in the census area as well.as the municipal boundaries. Figure 35 illust- 
rates the major trip desires in the Kitchener census area and shows fairly 
clearly the major commuter-shed sub-regions. Some of the linkages illust- 
rated in Figure 35 are also due to the timing of development. 


Windsor 


Figure 36 identifies three clusters of census tracts in the Windsor census 
area. The CTs in cluster II to the east have strong linkages with CBD 
employment and to a lesser extent with the automobile manufacturing plant 
(CT 25). The CTs within cluster I are more strongly linked to employment 
in CT 25. The dendrogram shows that cluster III consists of 13 CTs only 
three of which are illustrated in Figure 36 with the remaining 10 being 
located in the townships outside of the city boundaries. The dendrogram 
illustrates that there is a significant increase in the error sum of squares 
when cluster I and II are merged and that there is a sharp increase when 
these are merged with cluster III. The principal determinants of trip 
linkages in Windsor are socio-economic factors, the timing of development 
and to a lesser extent the magnitude of employment in CT 25. Figure 37 
illustrates the major trip desires in Windsor. 


London 


Figure 38 illustrates four clusters of census tracts for London. Cluster I 
embraces the CTs on the northeast fringe of the area which are newer re- 
sidential areas with commuting linkages to the CBD (CT 22). Cluster II 
consists of the CTs to the east of the CBD which also have strong linkages 
with the CBD and to employment opportunities in CT 27. Cluster III embraces 
CTs to the west of the CBD and south of the river. These CTs also have a 
strong orientation to jobs in the CBD. Cluster IV is to the north of the 
river and these CTs tend to have strong linkages with employment at the 
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FIGURE 38. Census Tract Clusters for London 
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University of Western Ontario and to the CBD. The dendrogram shows that 
clusters I and II merge and clusters III and IV merge and that there is 
not a sharp increase in the error sum of Squares when these merges occur. 
Figure 39 illustrates the major trip desires in London and demonstrates 
that commuting patterns are strongly dominated by the CBD employment and 
to a much lesser extent by the employment at the University of Western 
Ontario. The rivers and socio-economic factors also have an important 
impact on commuting patterns. 


>t. Catharines 


Figure 40 illustrates the three major clusters of census tracts for the 

St. Catharines census area and this diagram shows that the cluster boundaries 
follow the municipal boundaries of St. Catharines, Niagara Falls and Welland. 
The dendrogram inthe appendix illustrates that the CTs within each municipality 
cluster at low error sums of squares but that a very large increase in the 
error sum occurs when the three clusters are merged. Commuting patterns 
within each of the municipalities are determined principally by socio- 

economic factors. Figure 41 shows the trip desires for the area and these 
reflect the self-containment already mentioned in the dendrogram. 


Hamilton 


Figure 42 illustrates the census tract composition of five clusters for the 
Hamilton census area where three are located within Hamilton. and two are 
within Burlington. The CTs in cluster I are located on the Niagara Escarp- 
ment and have strong linkages with the heavy industrial areas. Cluster II 
consists of CTs in the west end of Hamilton which have higher socio-economic 
Status and these CTs have strong linkages with CBD employment and employment 
at McMaster Universtiy. Cluster II in the east end consists of lower status 
CTs which have strong linkages with the steel mill areas. The dendrogram in 
the appendix illustrates that a significant increase in the error sum of 
Squares occurs when clusters I and II are merged with cluster III and that 
very large increases occur when the Hamilton clusters are merged with the 
Burlington clusters. In Burlington the CTs in cluster IV have strong link- 
ages with Hamilton while those in cluster V are more strongly linked with 
employment opportunities in Burlington. Figure 43 illustrates the major 
trip desires in the Hamilton census area. Job concentrations, topographic 
features and socio-economic factors are the primary determinants of com- 
muting patterns in the Hamilton census area. 


Ottawa 


Figure 44 illustrates eight census tract clusters for the Ottawa census 
area six of which are within the Ottawa-Carleton region and two of which 
are in Hull. Clusters II and IV embrace the older residential areas of 
Ottawa while the remaining six on the southside of the river consist pri- 
marily of the newer residential areas. The dendrogram inthe appendix shows 
that clusters I and II which consist of CTs in a southerly corridor cluster, 
clusters III and IV in the east end cluster and clusters V and VI in the 
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west end cluster. The next level of clustering illustrated by the dendrogram 
is the integration of clusters I, II, III and IV in the eastern half of the 
areas. There is a significant increase in the error sum of squares when all 
clusters on the southside of the river are joined and a sharp increase when 
clusters on both sides of the Ottawa River are joined. Commuting to the 
central area tends to dominate the commuting pattern as shown in Figure 45. 
In addition the timing of development has a strong influence as do socio- 
economic characteristics and topographic features. 


2.4 Some Generalizations About Commuting Patterns 


The previous paragraphs have outlined the census tract composition 
of the major clusters of residential areas that have similar destination 
patterns. The principal determinants of these commuting patterns have been 
highlighted for each census area under each of the five broad categories 
identified in Section 2.3. Table 4 summarizes the major determinants of 
commuting in each of the fifteen census areas in terms of each of these 


areas. 


Seven of the census areas as defined by Statistics Canada embrace 
two or more separate municipalities. In almost all cases major cluster 
boundaries were coincident with the municipal boundaries indicating that 
commuting patterns tend to be self-contained within municipalities. Thunder 
Bay, Kitchener, St. Catharines and Hamilton provide very good examples of 
this type of effect. 

Topographic features such as rivers also have an important in- 
fluence on commuting patterns. The best examples are provided by London, 
Hamilton and Ottawa. Ottawa is a special case in that the Ottawa River is 
also a provincial boundary. 

One of the most important influences on commuting patterns is 
the timing of development where residential areas and employment zones 


that expanded during the same time period tend to have strong linkages. 
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The best example is provided by Guelph where the residential areas and 
employment zones that grew rapidly during the 1966-1971 period were all 
located on the periphery. The commuting patterns that emerged as a con- 
sequence of this phase of development have a distinctly different character 
to the centrally focussed patterns that developed earlier. 

Perhaps the strongest influence on the commuting structure of a 
community are the socio-economic factors that influence dwelling unit 
location decisions and employment location possibilities. The importance 
of socio-economic factors is analyzed in more detail in subsequent chapters 
of this report. 

The final factor identified in Table 4 is the impact that specific 
employment concentrations have on commuting patterns. Communities which 
have large concentrations of employment in one area either because of one 
industry or a collection of employers in that area develop particular types 
of commuting patterns that must be recognized in any modelling attempts. 
Sarnia, Sault Ste. Marie, Oshawa, Sudbury, Windsor and Hamilton provide 
examples of areas with strong concentrations of industries in particular 
areas and the commuting patterns to these areas tend to dominate the home 
to work linkages patterns. London and Ottawa are examples in which employ- 
ment in the tertiary sector pate tee the economic base of the community 
and this employment type tends to be concentrated in the CBD. 

The five determinants identified above and the associated clusters 
described in the previous sections provide a basis for identifying calibration 
sub-regions within each of the census areas. The particular sub-regions used 
in the calibration of multi-parameter gravity models are described later in 


this report. 


Some understanding of the heterogeneity in commuting patterns in 
each of the census areas is provided by Table 5 ih which the maximum value 
of error sum of squares is reported for each census area along with the 
number of spatial linkages in each area. Clearly this magnitude will in- 
crease with increasing size of urban area but deviations from the general 
trend may be clearly identified. For example in the smaller urban areas 
the higher than average maximum error sum of squares for Sarnia, Sault Ste. 
Marie and Kingston may be noted. In the next size group Thunder Bay, Oshawa 
and Kitchener provide examples of census areas in which commuting hetero- 
geneity leads to sharp increases in the error sums of squares in the final 
stages of the clustering process. Finally in the three largest census 
areas analyzed the error sums of squares also increased sharply in the 
final clustering stages primarily because residential zones from separate 


municipalities were being forced together. 


TABLE 5. 
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Ontario Census Areas 


Census Area 


Guelph 
Peterborough 
Sarnia 
Brantford 
Sault Ste. Marie 
Kingston 
Thunder Bay 
Oshawa 

Sudbury 
Kitchener 
Windsor 
London 

St. Catharines 
Hamilton 


Ottawa 


Number of Home- 
Work Linkages 


19,990 
20, 400 
23,760 
24,195 
24,740 
26,770 
33,130 
33, 740 
45.580 
81,095 
74,455 
97,195 
91,815 
152,265 


201,015 


Maximum Value of Error Sum of Squares for 


Error Sum of 
Squares 


SIRE 


a 


APPENDIX TO CHAPTER 2 


DENDROGRAMS FOR ONTARIO CENSUS AREAS 


The dendrograms for the fifteen Ontario census areas 
are included in this appendix. These dendrograms show the 
the hierarchical structure of residential census tracts which 
have been grouped on the basis of similarity in destination 
zone characteristics. 

The dendrograms are all plotted to the same error 
sum of squares scale along the ordinate in order to illustrate 
the degree of heterogeneity in commuting patterns. Error sums 
of squares are not directly comparable since the magnitude does 
increase with the number of census tracts used. The official 
census tracts only are included in the maps illustrated in 


this appendix. All census tracts are included in the dendrograms. 
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CHAPTER 3 


GRAVITY MODEL CALIBRATION AND 
GOODNESS OF FIT STATISTICS 


The calibration procedure used to aattnere the parameters of 
alternative forms of the gravity model is described in this chapter. In 
addition a variety of goodness of fit statistics which may be used to 
assess the qualities of any calibrated gravity model are described in 
the second part of the chapter. The relative effectiveness of these 
alternative statistics are also examined for a range of Ontario census 


areas. 


3.1 Basic Gravity Model Structure 


The trip distribution model which is calibrated in subsequent 
chapters to the 1971 journey to work data for all census areas is a 
doubly constrained version of the gravity model with a negative exponential 


deterrence function of the following form: 


* 
T., = A,B,0.D, exp(-8_c,,) (3) 
14 5 ae | ae Ba | s ij 
* e 2 e 
where T.. = the model estimated trip interchanges between zones 1 
1j e 
and is) 
0, = trip productions of zone i 
ae = trip attractions of zone j 
8 = travel deterrence function parameter specific to all the 
5 origin zones in a given calibration sub-region s 
or = the travel "costs" (distance, time or some combination) 
J 


for a trip between zones i and j 


mee Lf 


and 
-1 

A ie B.D. exp (—B c,.,) (4 
i UL yMeewers a ince ) 
B. = [Z A,O, exp(-8 c¢ ‘ike (5) 
4 i eS 2 Sane 

two balancing factors which ensure that the trip end constraints 
Sheed ER hagas (6) 
i ae | 

J 

and 
Doman Te. (7) 
a i J 


are satisfied. 


3.2 Model Calibration Technique 


The magnitudes of the balancing factors A, and ae are estimated 
by assuming initially that all S are equal to 1 and calculating all of the 
A, for a specific set of B values. These A, magnitudes are then sub- 
stituted into equation (5) to calculate a new set of os magnitudes. This 
process continues in an iterative manner until the product magnitude nig 
converges for all ij pairs. Convergence is assumed when the difference 
between the products of rhe for two successive iterations is less than 
1.0 x 10°’ 

The balancing factors A; and S may be set to 1.0 in order to 
produce, respectively, the production-constrained and attraction-constrained 
versions of the gravity model. These forms of the gravity model are cali- 
brated for selected census areas. 

The gravity model is calibrated by establishing the magnitudes of 


B which minimize the sum of the absolute differences between the observed 


and model-estimated ordinates of the trip length frequency distributions. 
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These frequency distributions are specified in 1 km class intervals and 
where more than one calibration sub-region is used, the observed and simu- 
lated frequency distributions are calculated separately for each B using 
only trips originating from those zones to which each parameter is specific. 
The model parameters are estimated using a golden section search 


technique [5]. With this procedure an initial estimate is made of the maximum 


and minimum limits for the parameter, Bo and Boo - Iwo new values are then 
chosen such that: 
She ge ChE Cre ae oe 
Pig OsBLohen ae Eqi2 ie -Poo O) 


The magnitudes 0.382 and 0.618 are derived from (T-1)l and 1/If where [ = 1.618 
the asymptotic interval reduction factor of the Fibonacci search method. The 
model estimated trip length frequency distributions are calculated using 


B 


ll and Bio and which ever magnitude gives the worst fit is chosen as the 


new limit of search on that side of the optimum. For example, if Bio were 


the worst the next two parameter magnitudes would be: 


pan SR Sap Rann 7 ten 
Ron 7 DOES (Bs 4 = Bat) eBay 
However, because of the properties of the Fibonacci numbers 8 = 6 and 


ide Ze 


only one new parameter and therefore one new frequency distribution need be 
calculated during each interation. 

Figure 46 illustrates the steps followed in this search procedure 
and convergence is obtained when the difference between two consecutive 


values of 8 is less than some arbitrarily set number and whichever of 


FIGURE 46. 
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Golden Section Search Procedure 


Ba, 83) 


B\2 Boo 


Bi2 


et? 


these g magnitudes is best according to the calibration procedure is 
taken as the parameter magnitude, 

The adjustment process for A, and a, and the estimation of 
Be are conducted in an inter-dependent way. The A, and Ba are ad- 
justed for each B before the calculation of the model estimated trip 
length frequency distribution. The logic of this process is illustrated 


in Figure 47. 


3.3 Alternative Model Statistics 

The calibration criterion used in the estimation of the gravity 
model parameters is the minimization of the differences between the 
observed and model estimated trip length frequency distributions. How- 
ever, the observed and model estimated trip matrices may be compared in 
many ways and the following sections identify some alternative statistics 


calculated for each of the gravity models calibrated. 


3.3.1 Mean Trip Length Statistics 


Three mean trip length statistics are reported for each of the 
calibrated models described in subsequent chapters of this report and 
these are the observed, simulated and expected mean trip lengths. The 
observed and siinilated mean trip lengths require no explanation and the 
expected mean trip length is calculated from the trip matrix estimated 
by equation (3) where 8 = 0. It is the mean trip length that would be 
expected to result from the spatial structure of a community if there 
were no travel deterrence effect. The expected trip length should be 
much higher than the simulated or observed where the difference magnitude 


would be related to the influence of the 8 magnitude. 
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FIGURE 47. Calibration Procedure for Doubly Constrained 
Gravity Model 
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Mean trip length statistics are also calculated for each census 
tract with one set relating to the trips originating from a zone and a 
second seCirelative to the trips terminating in a zone. The census tract- 
specific statistics allow census tracts with inadequate performance to be 
identified as well as assisting in the assessment of the effectiveness of 


both trip end balancing and multiple parameter modifications. 


3.3.2 Chi-Squared Statistic 
The chi-squared statistic is calculated from the model output in 


the following way: 


2 Oo Tah 
ij hg 
1j 
where T,., = the observed trip interchange magnitude between zones i 


and j 
The chi-squared statistic is also calculated for the individual origin and 


destination zones. 


Fe 
ae , 2 Se tale sb 
Origin-Specifie y» = 2 Ss a se Crip) 
i Th 
Chee Taps 
: : wre 2 1j 1j 
Destination-Specific y° = f z (12) 
j T 


This is a very sensitive statistic for use in comparing the 
simulated and observed trip matrices but there are some difficulties with 
its use. Interchange magnitudes of less than five in either matrix violate 
the assumptions underlying the chi-squared statistic. In addition the 


generally poor fit of the model in a statistical sense prevents this 
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Statistic from being used in conventional hypothesis testing. 


3.3.3 Phi-Statistic 
The phi-statistic is used in information theory and is calculated 
from k 


g==22 Tepito as (13) 
6 


The larger the absolute value of phi the poorer is the degree of fit of the 
model. Since each term in the calculation is multiplied by the observed we 
magnitude this statistic is most sensitive to errors in large interchanges. 


This statistic may also be calculated for the individual origin and destin- 


ation zones. 


3.3.4 R-Squared Statistics 


Two R-squared statistics may be calculated in the following ways: 


: Divan Gierea 9a 1s 
R nia NE NEN aa (14) 
Doe CTyayen Ta.) 
, i ac 
Int 
and 
x 
9 Fag ® (Ti. m= T, 4) 
|e ae ced sag A ia ens ay (15) 
adj = 2 
DY OreGl ey eal ae) 
ee ij 17 
O.D 
h Gt aera 
where ij T 


the expected number of trips between a zonal pair if there 
were no effect of travel deterrence 


Equation (14) compares the differences between the observed and 


estimated trip interchange residuals with the differences between the observed 


—nLO6 T= 


. F : ‘ 2 
and mean trip interchange magnitudes. In this equation R° measures the 


# 
» B, and 8 on the degree to which T,, is an 
1 j 1j 


improvement on lhe - On the other hand equation (15) tends to measure 


effects of O Dd. > A 


i b 
* 

the direct effects of 8 , or the degree to which re is an improvement 

on ti - Because of this the Ee magnitudes will tend to be lower than - 


the R magnitudes, 


3.3.5 Likelihood Statistics 


Three likelihood statistics may be calculated and these are: 


lee 

Lops TP : T,, log, = (16) 
ie 

Le-o = : yi a log. a ly) 
* 
a 

Loam = : Hi 114 log, “er (18) 


Equation (16) calculates the highest likelihood magnitude that can exist for 
a given urban area since it operates on the observed te magnitudes. 
Equation (17) calculates the magnitude of the likelihood that would be 

given by a trip distribution model without interaction deterrence effects. 
Finally, equation (18) calculates the likelihood foe a model with optimum 


8 and this likelihood magnitude should fall between the two extremes. 


3.3.6 Other Error Statistics 
A variety of other error statistics may be calculated and these 
include the following traditional measures. 
Mean Error 
MERR = ££ (T 
ij 


* vy 
ae T,,)/n (19) 


a LOT t= 


Mean Absolute Percentage Error 


* 
ie ; 
MABSPC = £ aan (20) 
a he | ch 


ij 


Standard Deviation of Residuals 


* 
(Ts, - Ty — MERR : 
SDRES = [oo yt = } pl2 (21) 
5 A n -l1 
Total Absolute Error 
* 
TABSERR = CE ET. = T), (22) 
eae i 
tej 
Mean Absolute Error 
MABSERR = TABSERR/n< (23) 


3.4 Detailed Analyses of Chi-Squared and Phi Statistics 


Previous analyses have shown that the chi-squared and phi statistics 
are very sensitive when used in the trip table context and large magnitudes 
are obtained. In order to examine further the behaviour of chi-squared and 
phi statistics Table 6 was generated in which the rows represent observed 
magnitudes of ve from 25 to 500 trips in increments of 25 trips and the 
columns represent simulated a magnitudes over the same range. The 
entries represent the contribution to chi-squared for each na versus ia 
comparison. The entries above the diagonal represent over-estimation residuals 
while the entries below the diagonal represent under-estimation residuals. 


Inspection of the entries in the matrix demonstrates clearly the asymmetric 


character of chi-square with the under-estimation errors being more important 
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than the over-estimation errors. 

Figure 48 shows the variation in the chi-square magnitude with 
ot4 Magnitude for constant error levels of +50, +100, +200. The diagram 
illustrates very clearly the rapid increase in the chi-square magnitude 
as the under-estimation error approaches the value of tes » that is as 
the value of vi approaches zero. This indicates the high sensitivity 
of chi-square to under-estimation errors when is is small which is a 
frequent occurrence in most simulated trip matrices. 

Table 7 shows the phi statistic magnitudes for different com- 
binations of observed and simulated trip interchange magnitudes. Inspection 
of Table 7 shows that the phi statistic is also asymmetric but less extreme 
than the chi-squared statistic. Figure 49 shows the variation in phi with 
the observed trip interchange magnitude for different error magnitudes. 
Although phi has a greater sensitivity to under-estimation than over-estin-— 
ation this is not nearly as severe as was evident in the analysis of chi- 
Squared. It should be noted also that the rapid increase of the statistic 
when ag approaches zero is most marked for large values of TG . The 
phi statistic is insensitive to errors involving small interchange values. 
Inspection of the over-estimation values reveals this property even more 
clearly. The phi values increase as oe increases despite the fact that 
the relative size of the error as compared to the observed value decreases. 
This is a very useful property in a statistic being used to evaluate trip 
distribution models. In most cases the observed matrix has been synthesized 
using a limited sample which often results in low levels or reliability for 
small interchange values. The phi statistic automatically accounts for some 
of this unreliability by discounting the effect of errors in small inter- 


changes. 
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FIGURE 48. Variation in Chi-Squared Magnitude with 
Observed Trip Interchange Magnitude 
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FIGURE 49. Variation in Phi Magnitude with Observed Trip 
Interchange Magnitude 
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The analysis of the response characteristics of the chi-squared 
and the phi statistics at the interchange level indicates that each re- 
sponds in a slightly different manner when the error changes. Both mea- 
Sures have a greater sensitivity to under-estimation errors than to over- 
estimation errors with this characteristic being less severe for phi. 

Phi is relatively more sensitive to over-estimation errors in large inter- 
changes with chi-squared more sensitive to allerrors in small interchanges. 
For under-estimation errors both statistics are more sensitive to an error 


of a specific magnitude on a smaller interchange. 


3.5 Goodness of Fit Statistics for Actual Census Areas 

The following sections of this chapter examine the behavioural 
of the goodness of fit statistics described in Section 3.3 for selected 
Ontario census areas. The behaviour of the statistics is examined both 
within individual census areas and between census areas. It is not pos- 
sible to approach this problem of the identification of the best goodness 
of fit statistics in the traditional statistical sense. The magnitudes 
of the statistics are influenced by the number of zones, the total number 
of trips and the characteristic probability distributions of the observed 
and estimated trip interchange magnitudes. A simulation approach has been 


taken in order to evaluate the performance of the alternative Statistics. 


3.6 The Simulation Method 

Eight Ontario census areas were selected for the detailed evalu- 
ation of model statistics and Table 8 lists these areas along with the 
population, number of census tracts and total number of observed trips. 
The areas selected represent a range of population sizes and spatial 


structures. 
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TABLE 8. Characteristics of the Eight Ontario Census Areas, 1971 
Census Area Population Number of Number Mean 
In Official Census of Interchange 
Census Tracts Tracts Trips Magnitude 
Oshawa EAS rs a 22 33, 740 70 
Thunder Bay 108,411 25 SSL 53 
Kitchener PLZ Os 45 81,095 40 
Windsor 211,494 56 74,455 24 
London 22.22 59 97,195 30 
St. Catharines 2202262 53 91,815 33 
Hamilton 437,554 109 152,265 LS 


Ottawa 506, 316 120 201,015 14 


AG Ge 


Simulated trip interchange matrices were calculated by multiplying 
the observed trip interchange entries by a random percentage error where these 
errors were generated from a rectangular distribution with mean of zero and 
a specified range. Six percentage ranges were used and these were ten, 
twenty-five, fifty, seventy-five, one hundred and one hundred and fifty 
percent. The random number generator GGU3 from the IMSL subroutine library 


was used and a simulated matrix calculated from the following equation: 


ieee eet Gel et RND. FACT Cy) 
1j 1j 1j 
where 6 = a random number taking values of +1 or -1 
RND = a random number where O < RND < 1.0 
FACT = allowable percentage error divided by 100 


The same starting value or "seed" was used for each set of simulations for 
each city and the estimated trip interchange matrix was simulated three times. 
The use of the same starting numbers made it possible to evaluate 
each statistic with respect to the variations in the allowable error only. 
The multiple simulations were used to ensure that these evaluations were 
robust and were not dependent upon a particular sequence of random numbers 
that happened to be generated. The simulated trip interchange matrices were 
then adjusted so that the row and column totals equalled the observed labour 
force and employment magnitudes. In this sense the simulated trip interchange 


matrices were similar to those produced by a doubly constrained gravity model. 


3.7 Behaviour of Alternative Statistics 


Figure 50 shows the variation in the Re and R magnitudes 


dj 


with allowable error percentage for the eight census areas. The tabulated 
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magnitudes are presented in the appendix of this chapter. This diagram shows 
that in all cases the eee Statistic is lower than R- and more sensitive 
to error increase. It is interesting to note that even when the allowable 
percentage error is 100 percent the Ro Magnitudes are greater than about 
0.7. The diagram also illustrates that R? is not a very sensitive statistic 
for comparing the observed and simulated trip matrices in the larger cities. 
This is also the case for the Bae in the larger cities except for Ottawa 
where it continues to decrease with increasing error percentage. 

The star marked on each plotted ne line indicates the R- magni- 
tudes obtained for the gravity models calibrated for each of the eight census 
areas. The locations of these points would suggest estimated trip matrices 
with errors in the range of 60-70 percent except for Windsor, Hamilton and 
Ottawa where error magnitudes of around 90 percent are indicated. 

Figure 51 shows the variation in the chi-squared statistic with 
percentage error for the eight census areas and the tabulated values are 
presented in the appendix. This diagram illustrates that the chi-squared 
magnitude increases sharply beyond an allowable percentage error magnitude 
of 75 percent. Beyond an error magnitude of 100 percent the increase in 
chi-squared tends to level off and this reflects the fact that negative 
trip interchange estimates are not allowed and the under-estimations are 
constrained to a maximum of 100 percent. The stars shown on the diagram 
identify the chi-squared magnitudes for the calibrated models and these 
correspond generally to error magnitudes of 75-100 percent. 

Figure 52 shows the variation in the magnitude of the phi statistic 
with increasing percentage error. It is important to note that upstLoverroys 


magnitudes of 100 percent the phi statistic increases almost linearly with 
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FIGURE 51. 
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increasing percentage error. The phi statistic alsé appears to be more stable 
than chi-squared in that the ranking of cities is consistent across all error 
magnitudes except for Kitchener and Windsor. The phi magnitudes for the 
doubly constrained models estimated for each census area are shown on the 
diagram suggesting interchange error magnitudes of the order of 80-90 percent. 

Figure 53 illustrates the variation in the likelihood statistics 
with increasing percentage error for the eight census areas and the tabulated 
magnitudes are presented in the appendix. The relationships illustrated show 
that the likelihood statistic behaves very similar to the chi-squared statistic 
with the magnitude increasing very sharply beyond error magnitude percentages 
of about 75 percent. 

Figure 54 shows the standard deviation of the residuals with in- 
creasing error magnitude. This diagram illustrates that the increase is 
roughly linear with the exception of London where the magnitude decreases 
with a one hundred and fifty percent allowable error. A similar behaviour 
may be observed in the Ro statistic for London and this behaviour occurs 
only in the first of the three simulations and appears to be a result of 
the particular random number set generated. It is also interesting to note 
that the standard deviation decreases with increasing census area size and 
this is because the average size of an interchange is much larger for the 
smaller cities. 

Figure 55 shows the variation in the mean absolute error of a 
trip interchange magnitude with increasing trip interchange error. The 
mean absolute error behaves in a very similar way to the standard deviation 


of the errors. 
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3.8 Comparative Behaviour of Statistics 


Three criteria may be identified for evaluating the quality of the 
alternative statistics and these are: (i) sensitivity to changes in error; 
(ii) similarity in magnitude across cities with different sizes; and (iii) 
consistency with other statistics. 

The graphs of each of the statistics as a function of the maximum 
percentage error have shown that phi, the standard deviation of the residuals 
and the mean absolute error are equally sensitive to change in the error 
across all magnitudes of the error percentage, The chi-squared and likeli- 
hood statistics increase exponentially with increasing error magnitude and 
do not provide good discrimation at low error Magnitudes. The Re statistic 
had very poor sensitivity to error magnitude changes remaining at high levels 
for large changes in the error magnitudes, The Bs statistic is a more 
sensitive statistic for some of the urban areas but was insensitive in some 
urban areas particularly London and St. Catharines. 

With the exception of the R and ay Statistics none of the 
Statistics may be used to compare the goodness of fit of models across all 
urban areas since the magnitudes of the statistics are sensitive to urban 
area size, or the number of trip interchanges entries being compared. How- 
ever, some transformations of the Statistics are suggested in the following 
sections to allow comparisons across urban areas. 

The third criterion suggested above is consistency. The chi- 
Squared statistic appears to be inconsistent for large error magnitudes 
Since the error functions cross-over at large percentage error magnitudes. 
Also in one of the simulations for R? and standard deviation of the residual 


measures are inconsistent for London where both magnitudes decreased when the 


error percentage increased from 100 to 150 percent. 
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Table 9 summarizes the ratings of each of the statistics with 
respect to the three criteria. This would suggest that the phi statistic 
and the mean absolute error statistics are perhaps the best statistics to 
use. However, their principal deficiency is that they could be used across 


urban areas to compare model behaviour. 


3.9 Comparisons Across Urban Areas 


Chi-squared and phi statistics have been calculated for all of 
the Ontario census areas as well as for some stratified matrices for 
Kitchener, Hamilton and Ottawa using a maximum random error range of 75 
percent. Figures 56 and 57 show the magnitudes of chi-squared and phi 
plotted against the number of home to work linkages in each area. These 
relationships indicate underlying linear trends suggesting that a linear 
transformation of either measure may be used as an indication of model 
performance across urban areas. Perhaps the most useful measure, then, 
is the phi statistic because of its constant sensitivity to error magnitude. 

Figure 58 shows R? and the phi statistic divided by the number 
of interchanges plotted against the number of interchanges when the esti- 
mated matrix has been generated using a 75 percent error magnitude. This 
diagram shows that the transformed phi statistic is quite consistent across 
all urban area sizes while there is a band of variation of R? F 

The behaviour of these two statistics for the calibrated gravity 
models is shown in Figure 59. This graph shows that the transformed phi 
Statistic increases with urban area size indicating that model performance 
is poorer with increasing urban area size. 

Another possibility would be to use the standard deviation of the 


errors and the mean absolute error statistics to characterize the goodness 
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of fit of trip distribution models. However, an evaluation of these 
statistics across urban areas indicated that they had certain undesir- 
able characteristics. The analyses presented above would suggest that 
the transformed phi statistic is the best goodness of fit measure to 
use. This statistic is used only as evaluative statistic in this report 
but its properties would suggest that it should be used as a calibration 


criterion as well. 
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APPENDIX TO CHAPTER 3 


Error Statistic Summaries 


TABLE 3.1. 
Census Area 10% 
Let) 98 a98 
Oshawa E20 99 ee 
Le 99 se fe. 
120 99 296 
Thunder Bay LO 99 Bee: 
130 oo 98 
1207320 799 
Kitchener 1seQu YS oo 
1.0) B20 =99 
LO ate) as: 
Windsor 1-0 Bie, 2 hS 
0 S050 so, 
L0' E86 10 
London 1.0) 8b0 130 
1e0 AO 41 30 
1.0 BUG e130 
St. Catharines 1.0 1.0 1.0 
is 0° 1.0; 1.0 
1350p 0 oe 
Hamilton LO So 9 
LOO t20 29 
LO 160 ao 
Ottawa 1,0 7060 3 
1.0 99 99 
Cell entries are R? and R 
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Performance of R2 


Percentage Error 
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100% 


Census Area 


Oshawa 


Thuncer Bay 


Kitchener 


Windsor 


London 


St. Catharines 


Hamilton 


Ottawa 


TABLE 352+ 
10% 25% 
105 635 
Jos: 630 
101 605 

oF 543 
105 602 
113 684 
Z0L" 1578 
ZO" F522 
263 1479 
317 1706 
S26 BF 
289 1624 
aad ~“ESFZ 
362 =1908 
354 1954 
302 1647 
281 1605 
ZhOr TS 
636 3392 
689 3547 
662 3472 
881 4570 
854 4450 
866 4443 
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Performance of Chi-Squared 


Percentage Error 


50% 


2835 
2804 
Lose 
2483 
2780 
3092 
7026 
6880 
6621 
7524 
7765 
7258 
8213 
8323 
8659 
7269 
7326 
6773 
14757 
15526 
T5272 
20097 
19374 
T9374 


75% 


8645 

8220 

8214 

7635 

8731 

9450 
21431 
ZEZLO 
19961 
23004 
23862 
22460 
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25328 
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21562 
22803 
20702 
44635 
47796 
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59242 
59092 
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43838 
37447 
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35 387 
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45057 
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255 346 
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270496 
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275:722 


150% 


91660 
94628 
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27907 
99965 
144161 
241069 
1O2S2 1 
143492 
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16 3603 
142766 
188908 
19 3624 
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381554 
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594061 
359621 
414098 
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1558 
1543 
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1344 
1440 
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4332 
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Performance of Phi Statistics 


Percentage Error 
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TABLE 3:..4. Performance of Likelihood Statistics 


Percentage Error 


Census Area 


102% 25% 50% 754% 100% 150% 

re 309 1308 3411 6689 1379 

Oshawa 54 307 1292 3258 6305 10 349 
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CHAPTER 4 


AGGREGATE TRIP DISTRIBUTION MODELS 


This chapter exainines the goodness of fit changes for a 
range of aggregate gravity sede types. The explanatory qualities 
of production constrained, attraction constrained and doubly con- 
strained versions of the gravity model using distance as the inter- 
zonal cost measure are first examined. Changes in the calibration 
qualities are then examined when travel times and a combination of 
travel time and travel distance are used respectively as the cost 
measures. A variety of stratified trip distribution models are 


examined in Chapter 5. 


4.1 Production Constrained Models 

The previous studies of the census journey to work data [1] 
have shown that the doubly constrained version of the gravity model 
was best in simulating the observed home to work linkages. Production 
constrained and ab peapeton constrained versions of the gravity model 
have been estimated for eight Ontario census areas in order to provide 
a benchmark against which some of the modifications to the gravity 
model described in Chapter 5 may be judged. 

The basic form of the production constrained gravity model 


estimated in this study is: 


re 
I 


-1 (25) 


It 


where A Uy dD, exp (-6C, 5) ] 


=~ 139— 


and the terms have been defined previously in Chapter 3. 

Table 10 summarizes the characteristics of the production 
constrained models estimated for eight Ontario census areas, This table 
shows that the observed and simulated mean trip lengths are very close to 
each other in spite of the fact that the calibration criterion was 
minimization of the sum of the absolute differences between the ordinates 
of the observed and simulated trip length frequency distributions. The 
B magnitudes range from 0.152 in Windsor to 0.237 in the Kitchener census 
area and tend to decrease in magnitude with increasing mean trip length. 
The transformed phi magnitudes shown in the right hand column of 
Table 10 suggest that the goodness of fit of the production constrained 


models is equivalent to a randomly introduced error of 75-100 percent. 


4.2 Attraction Constrained Models 
ase SE eltton wonstrained Models 
The basic form of the attraction constrained model estimated 


injithis study is 


ie Baan, apc. 8 26 
ij eto jo ae Cig (26) 


ait 


where B Ly 0; exp (-8C, ,)] 


J 

The properties of the attraction constrained gravity models for 
eight census areas are summarized in Table 11. The £8 magnitudes are 
all smaller than the 8 magnitudes of the production constrained model. 
The 8 magnitudes tend to decrease with increasing mean trip length but 
at a faster rate than for the production constrained models. The phi 
Statistic and the transformed phi indicate that the goodness of fit of 


the attraction constrained models is not as high as for the production 


constrained models. 
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TABLE 10. Properties of Production Constrained Gravity Models 
for Eight Census Areas 


Mean Trip Length 


Census Area B (kms ) R2 Phi eRe 
Obs. Sim. 
Oshawa O,is7 4.10). 4424 028455 01259637 10938 
Thunder Bay SRE whe cle, wha bo Oe8s on Gl f B42 BSS 
Kitchener OZ 41 4.80 4.70 0.85 ISH/46 0.42 
Windsor 0.152 (Resk) Ps) 0.75 36,900 0.49 
St.Catharines 0.210 G20 OLE 0.94 41,801 0.46 
London ees Gods? 9 Gaile 0 2947 oe 5809 HO.0156 
Hamilton 0.160 S054 prod 0.76 mq BY 29 0857 


Ottawa O.202 7.63 7.48 0.74 130,;2648f G65 
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TABLE? 11. Properties of Attraction Constrained Gravity Models 
for Eight Census Areas 


Mean Trip Length 


Census Area B (kms ) RZ Phi = 
Obs. Sim. 

Oshawa Os133 4.10 4.22 0.83 Veh i480) 0.40 

Thunder Bay OF172 S337 te 5 O84 1 eeOD a te SO 

Kitchener 0.230 4.80 4.69 0.86 34,439 0.42 

Windsor 0.084 io Meee of eae O62 43,829)" 0.2.59 


St..Catharines Oe AO “6.50 6.06 0.94 44,385 0.48 
London 0.146 On er OU mes (9 2 5. 35 Omi) 
Hamilton Om? 8.05 Dod, Oc 07 sl03. Jom 0.66 


Ottawa One 2 7503 ie ee OO Oe UD te ye ys) 
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4.3 Doubly Constrained Models 

The doubly constrained gravity model estimated in this study has 
been defined in equation (3). The model has been estimated for the 
fifteen Ontario census areas and the properties of the calibrated models 
are summarized in Table 12. The observed and simulated mean trip 
lengths are close to each other for most census areas. The phi and 
transformed phi statistics shown in Table 12 indicate that the doubly 
constrained gravity model is superior to both the production constrained 
and attraction constrained versions of the model. While the 8 parameter 
tends to decrease with increasing mean trip length there is a great deal 
of variability in the relationship. The 8 parameter magnitudes for 
the doubly constrained model tend to be larger than the parameter 


magnitudes for the production constrained models. 


4.4 Travel Time and Time-Distance LULD Cost Functions 


Doubly constrained models have been calibrated for the Kitchener, 
Hamilton and Ottawa census areas using network travel times obtained from 
transport studies performed in these communities. The properties of these 
time based models are summarized in Table 13. This table indicates that 
the observed and simulated mean trip lengths are close to each other with 
the largest dipeerence being for the Kitchener CMA. The phi statistics 
indicate that the calibrated models are marginally inferior to those 
calibrated using network distances. It should be noted however that 
the difference decreases as the city size increases. This would indicate 
that distance-based networks are more effective for smaller cities and 
time-based networks more effective for larger cities with the cross- 


over point being in the range of 600,000 to 1,000,000 people. 


TABLE 12. 
Functions 

Census Area B 
Guelph 0.165 
Peterborough 0.288 
Brantford 0.206 
Sarnia U2 09) 
Kingston is 20D 
Sault Ste.Marie 0.131 
Oshawa Or 213 
Thunder Bay Orzo 
Sudbury 0.120 
Kitchener O29 
St.Catharines 0.214 
Windsor OSL 
London 0.180 
Hamilton O59 
Ottawa Ns 
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(kms ) 
Obs. 


Sh The 
35 JU 


4.20 


Mean Trip Length 


Sim. 


Phi 


5,742 
5,709 
7,114 
6,485 
8,281 
6,566 
11,008 
10,775 
20,876 
29,924 
37, 359 
34,625 
44, 360 
83,017 
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Properties of Doubly Constrained Gravity Models 
for Fifteen Census Areas Using Distance Cost 


TABLE 113. 


Census 


Area 


Kitchener 


Hamilton 


Ottawa 
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Properties of Doubly Constrained Gravity 
Models for Three Census Areas Using 
Travel Time Cost Functions 


Mean Trip Length 


8 (mins ) RZ Phi 
Obs. Sim, 
0.150 10.23 9.79 0.74 36,212 
0.121 13.38 13.48 0.76 88,430 


O7179 heyy! Ay 07/8) 121,352 


0.58 


0.60 


ase 


The doubly constrained model was also estimated using travel 


costs which represented various combinations of inter-zonal distance and 


travel time for the Kitchener census area. 


functions used 


ky 


tO 


IV 


where 


Cis 
1j 


a la yi 


time 


Cae 
ij 
(ae 
ij 


FACT 


FACT 


are: 


ars 


Poteeence ij time ar 


Chee hAGIn ob 


Pate ene: sy ne Tater, 


c, ,/Mean Trip Distance + tt, ,/Mean Trip Time 


[c, ,/Mean Trip Distance].FACT + [tt, ,/Mean Trip 


Time |.(1 —, FACT) 
"generalized" trip cost between zones i and j 


magnitude estimated for the doubly constrained 
model using network distances 


magnitude estimated for the doubly constrained 
model using network travel time 


network travel distance between zones i and j 


network travel times between zones i and j 


is an index used to weight distance and speed 
differentially depending on the average travel 
speed between zones i and j 


1 if average speed > 80 km/hr 


average speed - 10 


a if 10 < speed < 80 km/hr 


O if average speed < 10 km/hr 


The four generalized cost 


(27) 


(28) 


(29) 


(30) 
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With the generalized cost function specified in equation (27) 
the deterrence effects of distance and time are equally weighted. In 
equation (28) when the average speed is greater than 80 km/hr the 
deterrence effect is assumed to be all distance and when the average 
speed is less than 10 km/hr the deterrence effect is assumed to be all 
travel time. In equation (29) the distance and speed components are 
effectively normalized because they are divided by the area-wide mean 
trip distances and mean trip length and equation (30) has the same 
effect as equation (28). 

Table 14 summarizes the properties of the doubly constrained 
gravity models using the four different cost functions. These parameter 
magnitudes cannot be compared directly since the generalized travel cost 
units vary between each model. The observed and simulated mean trip 
lengths are close to each other but there is little change in the 
overall goodness of fit of the doubly constrained model. The models are 
of better quality than those using travel times but are marginally 
inferior to that estimated for the Kitchener census area using network 


travel distances as the travel costs. 


4. Comparisons.of Residuals Across Model Types 


It is also useful to examine the trip interchange residuals in 
some detail in addition to the goodness of fit statistics discussed in 
the previous sections. The trip interchange residuals for the Ottawa 
census area are examined for four separate model types in this section. 
The spatial distributions of residuals observed for Ottawa are broadly 
representative of the residuals for all census areas. 

Figures 60 and 61 show the over-estimation (simulated > 


observed) and under-estimation (simulated < observed) residuals for the 


TABLE 14. 


Cost Function 
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Properties of the Doubly Constrained 
Gravity Model for the Kitchener Census 
Area Using Different Cost Functions 


Mean Trip Length 


B R2 Phi 
Obs. Sim. 
0.490 hee oh 272 0.85 3iy123 
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production constrained gravity model estimated for Ottawa. Trips from 
the inner suburbs of Ottawa and the residential areas in Hull to the 
Ottawa central area employment zones tend to be over-estimated by the 
production constrained model. This tendency to over-estimate trips to 
the central area is directly related to the heavy concentrations of 
employment in the Ottawa CBD. The under-estimation residuals are 
primarily associated with trips to suburban employment locations in both 
Ottawa and Hull although there are some under-estimation residuals 
associated with the longer trips to the CBD. This latter type of 
residual is probably associated with the timing of development of certain 
residential sub-divisions and the growth of Federal Government job 
opportunities in the central area. 

Figures 62 and 63 illustrate the over-estimation and under- 
estimation residuals for the attraction constrained gravity model 
estimated for Ottawa. The over-estimation residuals tend to be 
clustered around the central employment area reflecting the tendency of 
the gravity model to allocate trips to the closest available opportunities. 
The under-estimation residuals are almost all associated with the long 
trips to the CBD employment opportunities and these residuals are a 
major cause of the poorer goodness of fit of the attraction constrained 
model. The frequency distribution of the labour force by census tract 
is much more uniform than employment and inspection of the 8-parameter 
magnitudes in Tables 10 and 11 shows that the 8-parameter magnitude for 
the attraction constrained model falls in order to compensate. 

Figures 64 and 65 show the over-estimation and under-estimation 
residuals Cora doubly constrained version of the gravity model. 

Figure 64 shows that the doubly constrained model tends to reduce the 


number of over-estimation residuals to the central area compared with 
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the production constrained model. On the other hand Figure 65 shows 
that the under-estimation residuals tend to increase for the doubly 
constrained model compared with the production constrained model. 

Figures 66 and 67 show the Over-estimation and under-esti- 
mation residuals for the doubly constrained version of the gravity model 
using network travel times rather than network travel distances. 

A comparison of the over-estimation residuals plotted in 
Figure 66 with those plotted in Figure 64 illustrates the tendency of 
the model to further over-estimate trips from some of the outer resi- 
dential areas of Hull. Comparisons of the under-estimation residuals 
in Figures 67 and 65 show that there is little change in the spatial 
pattern of residuals. 

The relative importance of the over-estimated and under- 
estimated residuals is further illustrated in Table 15 where the phi- 
Statistics are summarized for the four gravity model types. This 
table demonstrates quite clearly that the major source of error is 
the under-estimation of trip interchange magnitudes and that the superior 
behaviour of the doubly constrained model is due to a reduction in the 
under-estimation residuals. The table also demonstrates that the intra- 
zonal under-estimation residuals are a significant proportion of the 
total. However, in reviewing this decomposition of the residuals it 
should be recalled that the phi-statistic is more sensitive to under- 


estimation errors than to Over-estimation errors. 


4.6 Comparison of Residuals for Alternative Cost Functions 
or atternative Cost Functions 


A more detailed analysis of the doubly constrained gravity 


models calibrated for the Kitchener census area using alternative travel 
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TABLE 15. 


mS IN 


Relative Importance of Over- and Under- 
Estimation Residuals for Four Model 
Types for Ottawa 


Total Phi Intra-Zonal Phi 


Gravity Model 


Type 


Production 
Constrained 


Attraction 
Constrained 


Doubly 
Constrained 


Doubly 
Constrained 
With Time 


Over-Est. Under-Est. Over-Est. Under-Est. 


Residuals Residuals Residuals Residuals 
23,634 106,650 426 ZES9F 
alg PE 129,435 0 30,055 
21,049 98,460 873 16,622 


20,814 100,537 ps8, Los 33 
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cost functions is shown in Table 16. This table shows that there is 
little change between the different travel cost function types with 


the distance and distance-time travel cost functions providing the 


best goodness of fit. 
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TABLE 16. Goodness of Fit Statistics for Doubly 
Constrained Gravity Models Using 
Different Travel Cost Functions for 
Kitchener Census Area 


Total Phi Intra-Zonal Phi 
Travel Cost 
Function Type 
Over-Est. Under-Est. Over-Est. Under-Est,. 
Residuals Residuals Residuals Residuals 
Distance 8,178 225875 485 3,008 
Time 8,701 235/20 1,260 3,868 
hi 8,247 23,007 909 SSL / 
II 8,646 23; 107 teh st P33 fs)5) 
EEL 8,556 21,745 bylo2 2,943 


IV 8,752 2/3520 Lag 37231 
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CHAPTER 5 


STRATIFIED TRIP DISTRIBUTION MODELS 


This chapter examines the qualities of a number of stratified 
gravity models. The first group of stratified models examined in this 
chapter are multi-parameter gravity models in which separate 8-parameter 
magnitudes are estimated for each origin zone in production constrained, 
attraction constrained and doubly constrained forms of the gravity model. 
The second group of stratified gravity models discussed have separate B- 
parameter magnitudes estimated for from four to eight calibration sub- 
regions where the census tracts in each of these calibration sub-regions 
have been identified by the clustering procedure described in Chapter 3. 
The final set of gravity models evaluated in this chapter are stratified 


by socio-economic group. 


5.1 Multi-Parameter Production Constrained Models 

Table 17 summarizes the characteristics of the multi-parameter 
production constrained gravity models estimated for eight census areas. 
With these models separate 8-parameter magnitudes were estimated for 
each residential zone. The table shows that in all cases the area-wide 
observed and simulated mean trip lengths are close to each other for all 
census areas with the multi-community census areas of Kitchener and 
St. Catharines exhibiting the largest differences between the observed 
and simulated mean trip lengths. A comparison of the goodness of fit 


statistics in Table 17 with those presented in Table 10 for the single 
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TABLE 17. Properties of the Multi-Parameter Production 
Constrained Gravity Models for Eight Census 
Areas 


Mean Trip 


Census Area 8 Range Pirie Ro Phi as 
Obs Sim 

Oshawa 003.= 6472 ~ 4110 14616 290750 10,827 0.32 

Thunder Bay 032 si 7314 Sedt — Sevia) teas 10,7735. 0433 

Kitchener O77 = 479 4.80 4.59 0.89 30,420 0.37 

Windsor -015 = .414 7459 1 go 24 eaeoo 33,903) 0246 


St. Catharines sk27 = «367 6.50 5291, Gag? 36,164 0.39 
London -030 - .497 6.57 6.42 0.95 49 ,043 0.50 
Hamilton -032 - .489 8.05 7292 0.87 80,669 0.53 


Ottawa UG) = 497 97.563 7733 " OciSpesensoes 0.65 
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parameter production constrained models shows that the multi-parameter 
models all exhibited superior goodness of fit characteristics. The 
largest improvements were in the Kitchener and St. Catharines census 
areas. 

Figure 68 shows the spatial distribution of origin-zone 
specific 8-parameters for the Kitchener census area. It should be 
recalled that a 8-parameter magnitude of 0.237 was obtained for the 
aggregate production constrained model while the range obtained for 
the multi-parameter model was from .077 to 0.479. Figure 68 illust- 
rates that most of the fringe residential zones not located close to 
employment areas have lower than average 8 magnitudes while those 
for the inner residential areas closer to the large employment zones 
have much larger values. The calibrated 8 magnitudes for each re- 
Sidential zone reflects not only the mean trip length characteristics 
of the zone but the proximity of the zone to large concentrations of 
trip attractions. 

Figure 69 shows the spatial distribution of origin-zone 
specific f-parameter magnitudes for the London census area. The 8 
magnitude for the aggregate production constrained model was 0.177 
and the range illustrated in Figure 69 is from 0.030 to 0.497. Once 
again the lowest magnitudes are on the urban fringe and the largest 


values in zones close to the CBD and the University of Western Ontario. 


3.2  Multi-Parameter Attraction Constrained Models 
ee et traction Constrained Models 

Table 18 summarizes the properties of the multi-parameter 
attraction constrained models estimated for eight Ontario census areas. 


A comparison of the goodness of fit statistics in Table 18 with those 


- Toh 


qouayojTy UT uoftjnqtTazq std 7A9}eMeIe{-§J 


"89 qanola 


19¢- I6I- 


GOd- 


6Iz- = © QZL- 
pe ee ate 
6l0- 9¢2- G6I- 
ee glee zlz-- &Sz- 
ES" 2a. 7. cop Bole 
SI: 
ISb- 
E91: 


= 1657 — 


-043 °034 


‘084 


‘049 


‘073 
135 eye -075 
-075 ho) He 139 J 
eels 176 °448:1I91 


FIGURE 69. 8-Parameter Distribution in London 
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TABLE 18. Properties of Multi-Parameter Attraction 
Constrained Gravity Models for Eight 
Census Areas 


Mean Trip 
Length | 


Census Area 8 Range (iaing ) R Phi rh 
Obs Sim 

Oshawa «O%3 r -414 4.10 4.04 0.73 18,004 0.53 
Thunder Bay Ps oP ht ame ee FA LIS ibe ar esha Mie / $1 60235 
Kitchener « L2Ggs), 2492 4.80 4.61 +0985 99556240541 
Windsor -045 - .497 4209 Jag 0.64 44,577 0.60 
St. Catharines -103 - .409 6.50 5 .69 0.96 41,108 0.45 
London 073 = .497 ce 6.49 0.89 56, /08 OF es ts: 
Hamilton 3073 <7,497 8.05 1.90. “Vose we he 85 (meer, 74 


Ottawa MTS = OF i:G5, 6,87 0 lG0) Tae tig 90.72 
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in Table 11 indicates that the multi-parameter attraction constrained 

model has poorer goodness of fit characteristics than the aggregate model. 
The only census areas in which improvements occur are Thunder Bay, Kitchener 
and St. Catharines, the census areas with strong multi-community spatial 


structures. 


5.3 Multi-Parameter Doubly Constrained Models 


Table 19 summarizes the properties of the multi-parameter doubly 
constrained models estimated for the eight Ontario census areas. A com 
parison of the gocdness of fit statistics of Table 19 with those of Table 
12 indicates that the multi-parameter doubly constrained model is superior 
to the single parameter doubly constrained model. In addition the doubly 
constrained multi-parameter model is superior to the other two multi-para- 


meter models. 


ed Sub-Region Specific Models 


Table 20 summarizes the properties of production constrained 

gravity models estimated with separate 8-parameters for a number of 
individual calibration sub-regions for eight census areas. The cluster 
analyses described in Chapter 2 formed the basis for the calibration sub- 
regions isolated for each census area. This table shows that the observed 
and simulated mean trip lengths are close together for all census areas. 
A comparison of the goodness of fit statistics in Table 20 with those in 
Tables 10, 12 and 17 show that the sub-region specific model is superior 
to the single parameter production constrained model and inferior to both 
the single parameter doubly constrained and multi-parameter production 


constrained models. 


TABLE 19. 


Census Area 


Oshawa 
Thunder Bay 
Kitchener 
Windsor 

St. Catharines 
London 
Hamilton 


Ottawa 


~ h6S, = 


Properties of Multi-Parameter Doubly 
Constrained Gravity Models for Eight 
Census Areas 
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TABLE 20. Properties of Sub-Region Specific Parameter 
Production Constrained Gravity Models of 
Eight Census Areas 


Mean Trip 


Census Area B spk. RZ Phi Phi 
No. Range Obs Sim a 
Oshawa 4 -046 - .369 4.10 ae 0.88 IT343451 0.34 
Thee Bay 3 een te 0 Spo vee s mlta we Ao) 0.87 T1529 3) F OS 
Kitchener 3 SLoZe—) s oeoee 4, 50 4.55 0.87 Sep lay ae! Ve wale ah 
Windsor 3 <U098e— 7414 (ee) Tc4409 0.275 36 ,589 0.49 


st. Catharines 3 OOt = 223 6.50 Duo ae) 0.395 G1 223 0.45 
London 4 OY 3. 204 Onover 0 s01lem 0.94 DOLD Opa ve: 
Hamilton 7 -O5Ge=—2 229 8.05 S-0ome 0762 S05 SlOge (0 07 


Ottawa 8 e196 40:7 7205 7,45 0574 125.0 1585. 0.62 
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Table 21 summarizes the results obtained from the calibration 
of doubly constrained models incorporating specific 8-parameters for the 
same calibration sub-regions used for the production constrained models 
summarized in Table 20. A comparison of the goodness of fit statistics 
of Table 21 with those in Tables 12 and 19 show that the doubly constrained 
sub-region specific models performed at about the same level or marginally 
better then the single parameter doubly constrained models but were in- 


ferior to the multi-parameter doubly constrained models. 


5.5 Models Stratified by Socio-Economic Group 


A number of gravity models have been estimated using the strat- 
ified journey to work tabulations obtained from Statistics Canada and 
described in reference [1]. The three two-way classifications obtained 


are: 


Automobile Ownership 


The “number of automobiles in a household" variable formed the basis for 

this two-way classification where the following groups were identified: 

Group 1 -— no automobile 

- members of labour force from households with 0 automobiles 

- members of labour force who are non-heads of household from a household 
with 1 automobile 

Group 2 — automobile available 

- heads of households from a 1 automobile household 

- members of labour force from households with 2+ automobiles 


Tenancy 


Group 1 - owner occuped dwelling unit 
Group 2 - rental dwelling unit 


Period of Residence 


Group 1 - individuals residing in current residence for more than 5 years 
Group 2 - individuals residing in current residence for five years or less 
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TABLE c2 1} Properties of Sub-Region Specific Parameter 
Doubly Constrained Gravity Models of Eight 


Census Areas 


Mean Trip 
Length 


Census Area B fae) R Phi 
No. Range Obs Sim 

Oshawa 4 pep. oe OATES AROS) 4.03 0.88 Pt, 051 
Thunder Bay 3 073 ELI De dy eet age O87, 10.9127 
Kitchener fe 249 apni 2 4.80 4.50 0.90 30, 368 
Windsor 5 oa al his asi 7259 f0% 50.572, 34,594 
St. Catharines 3 ale: 0d MO eS0 wed. Holy 70,396 38, 384 
London 4 mil Ese ofS oe O55 7. £6.48 44 0294 44,273 
Hamilton 7 Ceo eld. 8.05 B16; -.0..183 835120 
Ottawa 8 SEG So ON) Jeb 3 F260 nee, O80 beet 11655139 
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Table 22 summarizes the characteristics of production constrained 
models estimated for members of the labour force without a car available 
and with a car available for the journey to work. The only goodness of fit 
Statistics reported in Table 22 are the coefficients of determination since 
these production constrained models were calibrated prior to the beginning 
of the research project described in this report. Inspection of Table 22 
shows that there is generally good agreement between the observed and simu- 
lated mean trip lengths for the two groups and that in most census areas 
the ®-parameter magnitude for the non-car owners is larger than for car 
owners. 

Table 23 summarizes the characteristics of doubly constrained 
models estimated for non-car owners in three census areas using both net- 
work distances and network times. Table 24 summarizes the information for 
the models calibrated for car owners. Since the zone systems and numbers 
of linkages are fewer for the stratified models it is appropriate to com- 
pare the transformed phi statistics in Tables 23 and 24 with those in 
Table 12. This comparison would suggest that the doubly constrained 
models stratified by car ownership status perform at a superior level to 
the doubly constrained single parameter models. The 8-parameter magni- 
tudes for the captive members of the labour force are significantly higher 
than for the non-captives for all three census areas reflecting the greater 
sensitivity of captives to travel costs. 

Tables 25, 26 and 27 summarize the characteristics of the 
stratified models calibrated for home owners and home renters. Tables 28, 
29 and 30 summarize the characteristics of the stratified models calibrated 
for periods of residence of greater than six years and for five years and 


less. The comparison of home owner and home renter behaviour indicates 
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TABLE .22. Production Constrained Gravity Models for 
No Car Available and Car Available Groups 
for Fifteen Census Areas 


Mean Trip Length 


2 
Census Area Group B (km) R 
Obs Sim 

No Car 0.10 3.40 3.60 Ons 
Guelph 

Car 0.14 3.90 3.90 Ont78 

No Car 0233 3.00 2.80 0.80 
Peterborough 

Can 0 i 3.40 2.24 0.88 

No Car 0.08 4.50 LO 0.68 
Sarnia 

Car 0.05 6.70 6.98 0.87 

No Car 8 ee es Saou 3.24 0.80 
Brantford 

Car ONzo 4.40 4.08 OFsv 

No Car Ue So50) 3562 0.79 
Sault Ste. Marie 

Car 0.20 4.90 4.98 0.92 

No Car ae lyg 4.30 4.56 saa 
Kingston 

Car OSLO 6.00 GLO 0.85 

No. Gare w07Z26 4.60 4.30 0.89 
Thunder Bay 

Car Onc 5.60 5.40 0.92 

No Car 0.15 3.90 4.10 D.74 
Oshawa 

Car Oe22 O10 4.10 0.80 

No Car eee 6.40 og he 0.72 
Sudbury 


Car 0.10 9.20 9.20 0.68 


TABLE 22. (continued) 


Census Area 


Kitchener 


Windsor 


London 


St. Catharines 


Hamilton 


Ottawa 


Group 


No Car 
Car 
No Car 


(ane 


No Car 


Gay 


No Car 


Car 


No Gar 


Car 


NouGat 


Gar 


Be ee 


Mean Trip Length 


Obs 


(km) 
Sim 


TABLE 23. 


Census Area 
and 
Deterrence 
Measure 


Kitchener 
Distance 


Time 


Hamilton 


Distance 


Time 


Ottawa 


Distance 


Time 


Properties of Doubly Constrained Gravity Models 
Calibrated Using Travel Time and Travel Distance 
For Three Urban Areas 


LoL 


22Of 


¥L69 


248 


Lb 


von 


Mean Trip 
Length 
(kms, mins) 
Obs Sim 


4.10 Seg. 


9.02 SETA 


Sayan! 5.80 


10.41 10.54 


6.15 Gets 


217,600 11.84 


: No Car Available 


Phi 


105239 


[UR EL 


“1,563 


23,577 


Do, 109 


357395 


TABLE 24. 


Census Area 
and 
Deterrence 
Measure 


Kitchener 


Distance 


Time 


Hamilton 


Distance 


Time 


Ottawa 


Distance 


Time 


~ 276 - 


Properties of Doubly Constrained Gravity Models 
Calibrated Using Travel Time and Travel Distance 
Car Available 


For Three Urban Areas : 


0.147 


8 a i 


0.158 


O5213 


Mean Trip 

Length 
(kms, mins) 
Obs Sim 
S08 49) 


10.61 10.48 


8.79 o2O0L 


13.70 14.13 


8. 34 8.51 


14.31 14.44 


Phi 


16, 310 


19,211 


So, Le 


42,207 


60,604 


37,245 


Cpe) 


0.46 


0.44 


TABLE: 25°. 


Census Area 


Guelph 


Peterborough 


Sarnia 


Brantford 


Sault Ste. Marie 


Kingston 


Thunder Bay 


Oshawa 


Sudbury 


Fa WI 


Production Constrained Gravity Models for 
Home Owner and Home Renter Groups for 
Fifteen Census Areas 


Mean Trip Length 


Group B (ian) 

Obs Sim 
Owners 0.14 3.60 3.80 
Renters OF 3.80 Seu 
Owners Oe23 SES) 3.10 
Renters 0.26 3.00 2.80 
Owners O07, 6.40 6.80 
Renters 0.02 5.70 5290 
Owners Ores) 4.20 Spele) 
Renters Oo2i 3.80 5.60 
Owners 0.16 4.60 4.60 
Renters 0.08 4.10 4.30 
Owners Ont3 5.80 Ds 00 
Renters 0.08 2220 awe) 
Owners VEPAG 5.40 on ae 
Renters Ox22 4.60 4.40 
Owners One 4.10 4.10 
Renters 0.16 3./0 Se 
Owners OSL2 8.60 8.40 


Renters 0.20 7.90 5.90 


TABLE 25. 


Census Area 


Kitchener 


Windsor 


London 


St. Catharines 


Hamilton 


Ottawa 


(continued) 


Group 


Owners 
Renters 
Owners 
Renters 
Owners 
Renters 
Owners 
Renters 
Owners 
Renters 
Owners 


Renters 


— GA = 


Mean Trip Length 


Obs 


4.70 


4.70 


7.10 


6.70 


6.60 


5.90 


6.30 


6.10 


8.30 


6.90 


8.40 


6.40 


(kin) 


Sim 


4.40 


4, 30 


6.30 


6.20 


5.30 


5.30 


eee, 


8.00 


6.70 


TABLE 26. 


Census Area 
and 
Deterrence 
Measure 


Kitchener 
Distance 


Time 


Hamilton 


Distance 


Time 


Ottawa 


Distance 


Time 


crip Fale haa 


Properties of Doubly Constrained Gravity Models 
Calibrated Using Travel Time and Travel Distance 


For Three Urban Areas 


0.170 


0.126 


Uietoa 


0.110 


Mean Trip 
Length 
(kms, mins) 


Obs 


err, 


HOR TS 


8.26 


ie 13 


8.43 


14.22 


Sim 


ec 2 


8529 


14.95 


: Home Owners 

R2 Pha 

0.89 ‘toga fae 3 
0.81 20,116 
0.91 Bye oe 
OFS3 44,016 
0.80 55,649 
04.82 50,916 


0.45 


TABLE 27. 


Census Area 
and 
Deterrence 
Measure 


Kitchener 
Distance 


Time 


Hamilton 


Distance 


Time 


Ottawa 


Distance 


Time 
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Properties of Doubly Constrained Gravity Models 
Calibrated Using Travel Time and Travel Distance 
Home Renters 


For Three Urban Areas : 


07192 


O75 125 


0.196 


65139 


Mean Trip 


Length 
(kms, mins) 
Obs Sim 
4.67 4.14 
9.90 9.97 
6.88 6.67 


L156 Leese 8. 


6.43 6209 


12.24 Ze 


Phi 


10,600 


11,5380 


21,039 


oe 500 


425722 


41,002 
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TABLE 28, 


Census Area 


Guelph 


Peterborough 


Sarnia 


Brantford 


Sault Ste. Marie 


Kingston 


Thunder Bay 


Oshawa 


Sudbury 


Production Constrained Gravity Models for 
6+ Years Resident and 0-5 Years Resident 


ah 6 er 


Groups for Fifteen Urban Areas 


Group 


6+ Yrs 
Geliy Vials: 
6FoYrs 
Oreoy Nias 
OF Yrs 
0-5 Yrs 
6+ Yrs 
0-5 Yrs 
6+ Yrs 
Y= S1.rS 
6+ Yrs 
= 3a. rs 
6+ Yrs 
Vaya te 
6+ Yrs 
OE Wares 
6+ Yrs 


0-5 Yrs 


Mean Trip Length 


Obs 


3.50 


3.20 


3.30 


6.00 


6.40 


4.10 


4.20 


4.50 


4.50 


5.40 


5.60 


Die G0 


220 


4.10 


3290 


(km) 


Sim 


TABLE 28. 


Census Area 


Kitchener 


Windsor 


London 


St. Catharines 


Hamilton 


Ottawa 


(continued) 


- 182 - 


Mean Trip Length 


Obs 


4.50 


4.80 


SEAN 


7.90 


6.10 


6.60 


6.00 


6.10 


7390 


7.90 


7.10 


7390 


(km) 


Sim 


ihe 
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TABLE 29, 


Census Area 
and 
Deterrence 
Measure 


Kitchener 
Distance 


Time 


Hamilton 


Distance 


Time 


Ottawa 


Distance 


Time 
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Properties of Doubly Constrained Gravity Models 
Calibrated Using Travel Time and Travel Distance 


For Three Urban Areas 
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TABLE 30. 


Census Area 
and 
Deterrence 
Measure 


Kitchener 


Distance 


Time 


Hamilton 


Distance 


Time 


Ottawa 


Distance 


Time 
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Properties of Doubly Constrained Gravity Models 
Calibrated Using Travel Time and Travel Distance 


For Three Urban Areas 


144 


151 


116 


- 164 


~119 


Mean Trip 

Length 
(kms, mins) 
Obs Sim 
4.83 4.68 


10.18 20.12 


vag B.S 


12559 13.08 


Veo 1 S.01 


13.84 ES. 85 


0.82 


: O-5 Year Residents 


Phi 


Lopcas 


1559637 


275905 


32,210 


54,848 


SE WF ie: 


0.41 


0.47 


0.44 
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that the differences are not systematic in terms of both the mean trip 
length comparisons and the 8-parameter comparisons. On the other hand 

the period of residence comparisons show that there are Significant dif- 
ferences between the two groups and that the 6-parameters in the three 
census areas for periods of residence of six years and greater are larger 
than for the shorter periods of residence. The transformed phi statistics 
show that this stratification has greater explanatory power than the 
single parameter model and is marginally inferior to the stratification 


based on car ownership status. 
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CHAPTER 6 


STRATIFIED TRIP GENERATION EQUATIONS 


Trip distribution models calibrated by socio-economic group 
have been described in the previous chapter. In addition reference [ay 
has described some preliminary trip generation analyses by socio-economic 
group for the Kitchener CMA. The special journey to work tabulations 
from Statistics Canada were not obtained for inclusion in reference [1] 
and this chapter describes the complete set of trip generation regression 
analyses conducted for all of the Ontario census areas. 

The first sections of the chapter summarize the trip generation 
and trip length characteristics of the socio-economic groups by residential 
tenancy and period of residence. The subsequent sections provide a more 
detailed analysis of the travel characteristics of non-car owners and 


car owners. 


6.1 Travel Characteristics by Tenancy and Period of Residence 


Table 31 summarizes the work trip linkage rates for the fifteen 
Ontario urban areas for dwelling unit owners and renters. The entries in 
this table are based on linkages to work only within the census areas arid 
do not include places of employment outside the census areas. 

The entries in this table show that the trip generation rate per 
household is higher for owners than renters with the exception of Kingston. 
The work trip linkage rate for owners varies from 1.04 in Kingston to bese h 
in Ottawa while the rate for renters varies from 0.90 in Thunder Bay to 


1.12 in Sudbury and in Ottawa. The higher average rates in owner occupied 


ns 


TABLE 31. Trip Generation Rates for Ontario Urban 
Areas Based on Tenancy 


Census Area Class Linkages Number of Trip Generation 
Produced Households Rate 
Guelph Own 14,790 11,590 Me Se: 
Rent 6,720 Ggb25 V.0% 
Peterborough Own 16,635 1277805 130 
Rent Sipe a) 6,005 0.93 
Sarnia Own 19,650 16,000 Bess 
Rent 6,675 ©" 6,635 POL 
Brantford Own 1/55 16,550 Lo 
Rent Pgh D 15300 So 
Sault Ste. Marie Own 20,235 175935 Los 
Rent 6,510 6,165 1.06 
Kingston Own Be 29 135220 104 
Rent 13,020 12,400 UR As 
Thunder Bay Own 30,585 Ae Sie As 1.30 
Rent ALIAS) 8,615 0.90 
Oshawa Own 25,965 oe ko Lia4 
Rent 105,320 1253 0.91 
Sudbury Own 30, 480 22,860 1.34 
Rent 33570 1635535 deo 
Kitchener Own 54,015 39,905 155.3 
Rent 30,105 26,645 1ST EAS 
Windsor Own Ga 3.55 52,240 a 
Rent 205245 215930 hes 
London Own 69,375 51,980 1353 
Rent 37,545 See) aw 1 be O17 
St. Catharines Own 77,805 63,895 be22 
Rent 22,905 25,050 Oa 
Hamilton Own Lr, 510 93,140 Lee 
Rent SME de. et 40 GL33 
Ottawa Own 120,780 85,595 4 


Rent 94,770 84,430 1,12 
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TABLE 32, Trip Generation Rates for Ontario Urban 
Areas Based on Period of Residence 


Census Area Class Linkages Number of Trip Generation 
Produced Households Rate 
Guelph 0-5 yrs. 2 ee Bo Bo 9,540 y ees HS, 
eo 10,410 8,665 20 
Peterborough 0-5 9, /695 8,650 VSD is 
6+ £25180 LO, 120 1520 
Sarnia 0-5 2 OU 10, 720 ua 
6+ 14,160 ids 9:35 tale) 
Brantford 0-5 12,000 10,685 reo 
6+ 14,910 £37230 1A? 
Sault Ste. Marie 0-5 11,430 dea 20 1.20 
6+ 5.300 11,570 Rep 
Kingston 0-5 "LF Lie 14,505 te ie 
6+ L370 10, 380 ee 
Thunder Bay 0-5 15,900 14,180 eee 
6+ 22,410 18,080 1324 
Oshawa 0-5 a7 g020 tf ALO 0.99 
6+ 195273 16,820 dsl 
Sudbury 0-5 26, 305 21; 075 io 
6+ 23574.) 18, 340 129 
Kitchener (y= 5 46,080 375999) Leal 
6+ 38,040 29,985 1.30 
Windsor 0-5 36 ,060 32,940 2209 
6+ 45,505 ALS145 Le 
London 0-5 SW et sie: bed 19S erat 
6+ 49,155 39 , 360 1325 
St. Catharines 0-5 43,065 SOLIS 17190 
6+ 57,660 49,765 dead6 
Hamilton 0-5 8262/2 73,943 pei: 
6+ 84,675 Peps20 Lior 
Ottawa 0-5 i226 00 93,595 24 


6+ 92,895 7iag425 1.30 
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dwelling units are due to a variety of factors which might include larger 
household sizes with greater numbers in the labour force. 

Table 32 summarizes the work trip generation rates for the 
fifteen urban areas based on period of residence categories. The gener- 
ation rates for the longer period of residence are uniformly higher than 
the 0-5 year category except for Kingston. The work trip linkage rate 
varies from 1.11 in Windsor to 1.32 in Sault Ste. Marie for periods of 
residence 6 years and greater and from 0.99 in Oshawa to 1.25 in Sudbury 
for periods of residence of 5 years and less. The lower trip rates for 
the shorter periods of residency reflect no doubt the smaller household 
Sizes of the recently established households. 

Table 33 shows the mean trip lengths of households classified 

by tenancy type and of households classified by period of occupancy. 
This table shows that the mean trip length to work of renters is signifi- 
cantly different from that of owners with the biggest difference being in 
Ottawa where many work trips from owner occupied dwelling units are from 
residences outside of the greenbelt. 

This table also shows that the average mean trip length to work 
from dwelling units that have penn occupied for five years or less is in 
most cases longer than the average trip from households that have been oc- 
cupied for six years or longer. This is reasonable since most of the new 
dwelling unit opportunities in urban areas are located on the periphery. 


The largest differences are in Windsor and Ottawa. 


6.2 Travel Characteristics by Car Ownership Status 


The most useful partitioned travel data available from the 1971 


Census for transport planning purposes are the Spatial interaction patterns 
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TABLE 33. Mean Trip Lengths for Socio-Economic Groups in Ontario 
Cities | 


Mean Trip Length 


(km) 
Occupancy Occupancy 
Owner Renter < 5 Years 6+ Years 
Guelph 3.6 3.8 a39 9 pv, 
Peterborough gus 30 oe | Si 
Sarnia 64 s) igi 6.4 6.0 
Brantford LS 3.8 4.2 aoa 
Sault Ste. Marie ge: ANS} 4.5 a5 
Kingston ee 52 5.6 Se 
Thunder Bay HG 4.6 52 be pres 
Oshawa ABS 057 3.9 ee 
Sudbury 8.6 729 Sey a ytd 
Kitchener 427 47 4.8 Ba 
Windsor a7 6.7 7.9 rep 
London 6-6 529 6.6 Gel 
St. Catharines aor 6.1 ‘eae f 6.0 
Hamilton 8.3 6.9 739 ‘eet 


Ottawa R22 6.4 V9 Pel 


2 A9Tb = 


by car ownership status. The ability to estimate transit captivity and 
non-captivity from dwelling unit information is a very useful capability. 

The regression analyses of work trip linkage formation rates 
reported in reference [1] demonstrated that the best prediction equations 
could be developed in terms of the dwelling unit composition of census 
tracts where the dwelling unit composition is expressed in terms of the 
number of single detached units and the number of attached dwelling units 
plus apartments. The use of the number of meee dwelling units as a 
third independent variable was found to yield Statistically non-significant 
partial regression coefficients for many of the census areas. An initial 
analysis of the trip productions by the two car ownership groups using 
three independent variables demonstrated similar behaviour. 

Tables 34 and 35 illustrate the prediction equations that have 
been developed to estimate the census tract amounts of captive and non- 
captive amounts of labour force, respectively, along with the statistical 
properties of these equations for the fifteen Ontario census areas. For 
the captive equations the partial regression coefficient of detached 
dwelling units varies from 0.156 in Kingston to 0.441 in Guelph and the 
t-statistics indicate that all of these coefficients are Significant at 
the 1 percent level. The partial regression coefficient of attached 
dwelling units varies from 0.331 in Oshawa to 0.639 in Guelph. Pooled 
regression equations are also shown for each census area size class and 
for all census areas combined. However the F ratios are all significant 
indicating that the pooled equations cannot be used instead of the individ- 
ual equations. Inspection of the two partial regression coefficients shows 


that labour force captivity generation from single attached dwelling units 


TABLE 34. 


Census Area 


1. Guelph 


2. Peterborough 


3h) soannia 


4. Brantford 


5. Sault Ste. Marie 


6. Kingston 


GROUP 1-6 
CHOW TEST 1-6 


7. Thunder Bay 


8. Oshawa 


9. Sudbury 


10. Kitchener 


GROUP 7-10 
CHOW TEST 7-10 


11. Windsor 


12. London 
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Census Zone Captivity Prediction Equations Using 


a Modified Dwelling Unit Composition 


=——y 


13 


16 


30 


-104 


+30 


10 


133 


-40 
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9G 
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200 
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O41 


OF 


apt 
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94 


37 
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nh hs, 
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F Ratio 


Ee wal 
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Ley gees 
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TABLE 34. (continued) 
b b 


Census Area 1 : 4 Re F Ratio DF 
Ce) GEES) 
ar 2 
13. St. Catharines 9 On7.8 0.537 0.93 
(6.32) (S205) 
14. Hamilton -38 0,205 0.536 Oo71 
(5.83) 7882 5') 
15.4. .0trawats 54 0.206 0.606 0.89 
(5.40) (18.91) 
GROUP 11-15 -13 0.219 0.589 0.86 
(12.80) (29.07) 
CHOW TEST, 11-15 i be-ws fete, eeu 
ALL CASES -5 OR224 0.586 0.88 
(19.82) (42.44) 
CHOW TEST — ALL CASES be 52. 42,298 
ZONE CAPTIVES =" a + by - Zone Single Detached Dwelling Units 


= i b, - Zone Single Attached Dwelling Units and 
Apartments 


a ae 


and apartments is from two to three times greater than for single detached 
dwelling units. | 

The non-captivity equations presented ih Table 35 show that the 
partial regression coefficient of detached dwelling units varies from 0.800 
in Oshawa to 1.267 in Ottawa while the cotati at of attached dwelling 
units varies from 0.185 in Brantford to 0.542 in Sudbury. The attached 
dwelling unit partial regression coefficients are not significant for 
Guelph, Thunder Bay, Windsor and St. Catharines. The non-captive co- 
efficients are two to three times larger for the detached dwelling units 
than the attached dwelling units. The table also illustrates that the 
pooled regression equations cannot be used in place of the equations for 
individual census areas. 

The labour force participation rates vary between the Ontario 
census areas and a second set of regression equations have been developed 
for the case in which the labour force participation rates have been 
factored to a constant magnitude across all census areas. This operation 
allows the relative importance of car ownership status by dwelling unit 
type to be compared more directly across census areas. The prediction 
equations developed in this way for both captives and non-captives are 
summarized in Tables 36 and 37, respectively, along with their statistical 
properties. Pooled regression equations for all census areas are also 
shown in the two tables. The F ratio for the pooled captive equation 
is not significant indicating that the pooled equation may be used to 
represent the fifteen census areas. The F ratio for the non-captive 


pooled regression equation is just significant at the 1 percent level. 


= Lo Se 


TABLE 35. Census Zone Non-Captivity Prediction Equations 
Using a Modified Dwelling Unit Composition 


Census Area by b, R F Ratio DF 
(ce 8. (en) 
1. Guelph 143 0.838 0.319 0.83 


2. Peterborough -12 0.989 0.262 0.94 
Be SSarnia -79 F073 02328 0.98 
4. Brantford 18 07953 Oet85 0.98 


3. Sault Ste. Marie 86 1.023 0.393, (0594 


6. Kingston BZ LALGS 0.247 0.98 
(7767) (4.04) 
GROUP 1-6 -22 1 O5G Oma02 0.95 
(S72) (7.09) 
CHOW TEST 1-6 3.364 A hs Wg 
7. Thunder Bay -165 1.068 -0.008 0.95 
(13.89) (-0.05) 
8. Oshawa 146 0.800 0.415 0.88 
(9.49) (3.35) 
9. Sudbury -94 14076 OF 542 0.91 
(12.58) ome) 
10. Kitchener SEE Loy Of 3538 0.86 
(10042) (2497) 
GROUP 7-10 -69 0.982 0.558 0.87 
(18.58) (8057) 
CHOW TEST 7-10 6.278 9,70 
11. Windsor 67 0.938 0.129 0.90 


12. London 130 1.039 0.320 0.89 


TABLE 35. (continued) 


Census Area 0, 
13. St. Catharines 118 
14. Hamilton -281 
15. Ottawa -25 
GROUP 11-15 =] 


CHOW TEST -11=95 


ALL CASES -20 


CHOW TEST - ALL CASES 


ZONE NON-CAPTIVES = a+ by 


+ Db 
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.020 


64) 


0.90 


F Ratio 


6.080 


3.7/9 


. Zone Single Detached Dwelling Units 


DF 


2 L5G 


42,298 


Zone Single Attached Dwelling Units and 
Apartments 


oe fe 


TABLE 36. Revised Captive Prediction Equations Using a 
Modified Dwelling Unit Composition 


b b 


C(M)A é * fs Ro 
t) C5 

Guelph -150 C2373 0.540 0.87 
(3.55) (5572) 

Peterborough LL D256 0.360 O25 .h 
(4.90) (5.06) 

Sarnia 14 0.180 0.408 OOF 
(77-10) Ciss12) 

Brantford 31 0.190 Oa 277. 0.91 
C5), 0.3) (9.41) 

Sault Ste. Marie -82 O22) 0.556 0.95 
(7.86) CLO) 

Kingston -26 0.134 0.532 D.95 
(5.32) (LEST) 

Thunder Bay 8 0.230 0.467 0.94 
(6.88) (7.19) 

Oshawa -124 OSL 0.309 OR87 
Gi 6) (5.16) 

Sudbury -32 Oe 216 0.443 0.86 
(S213) (9.16) 

Kitchener 40 OsbS3 OF534 0.78 
(2.80) C7 262) 

Windsor -89 0.244 05477 0.85 
(Oro) (9.95) 

London 14 0.208 0.424 0.94 
(8.18) (16.04) 

St. Catharines 8 QaL57 0.474 0.93 
(6.52) (8.15) 

Hamilton -33 OF 232 0.469 OekL 
(2030 (725) 

Ottawa 42 0.163 0.478 0.89 
(5-40) (18.91) 

ALL CASES zs 0.200 0.472 0.88 
) C27 746) (41.45) 
CHOW TEST F RATIO = 1.008 DF = 42,298 


Zone Captives/Aggregate Trip Production Rate = REVISED ZONE CAPTIVES 


REVISED ZONE CAPTIVES = a + by - Zone Single Detached Dwelling Units 


+ b, - Zone Single Attached Dwelling Units 
and Apartments 
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TABLE 37. Revised Non-Captive Prediction Equations Using a 
Modified Dwelling Unit Composition 


b b 
mye : te.) te) Re 
ty 2 

Guelph 121 0.709 0.270 0.83 
(O75) (1.89) 

Peterborough -10 0.846 0.224 0.94 
(133/22) (2565) 

Sarnia -68 0.923 Os2e2 0.98 
(21.96) (325)) 

Brantford 16 0.847 0.164 0.98 
(Z2e15) (2.89) 

Sault Ste. Marie 68 0.807 Ow2738 0.94 
(13.16) (2.45) 

Kingston 27 1,025 ORzi2 0.98 
(27587) (4.04) 

Thunder Bay -138 0.899 -0.007 0.95 
(13°89) (-0.05) 

Oshawa 137 0.749 0.388 0.88 
(9.49) Geis eh) 

Sudbury -74 0.848 Ov4a27 0.91 
(2258) (S590) 

Kitchener 22 0.943 0.289 0.86 
(10.42) (22943 

Windsor 61 0.852 Omid 7, 0.90 
CUse37) (15659 

London 106 0.847 0.260 0.89 
(12.69) (3.74) 

St. Catharines 104 0.902 OF0G2 0.96 
aL ds) (OL27) 

Hamilton -246 0.917 0.344 0.92 
(19.19) (4543) 

Ottawa -4] 0.999 0.288 0.86 
(16.06) pers 

ALL CASES -22 0.882 0.301 0.92 
(57,08) (1.5794) 

CHOW TEST F RATIO = 1.685 DF = 42,298 


Zone Non-Captives/Aggregate Trip Production Rate = REVISED ZONE NON- 
CAPTIVES 


REVISED ZONE NON-CAPTIVES = a + by . Zone Single Detached Dwelling Units 


+b, . Zone Single Attached Dwelling Units 
and Apartments 


Inspection of the captive partial regression coefficients of 
detached dwelling units shows that they vary from 0.134 in Kingston to 0.373 
in Guelph and those of attached dwelling units from 0.309 in Oshawa to 
0.556 in Sault Ste. Marie. The partial regression coefficients for the 
pooled equation indicate that for average labour force participation rates 
the captive generation rate'‘is 2.3 larger for attached dwelling units than 
detached dwelling units. 

Inspection of the non-captive partial regression coefficients 
of detached dwelling units shows that they vary from 0.709 in Guelph to 
1.015 in Kingston and for attached dwelling units from 0.104 in Brantford 
to 0.427 in Sudbury with the coefficients for Guelph, Thunder Bay, Windsor 
and St. Catharines being non-significant. The pooled equation suggests 
that the non-captive coefficient of detached dwelling units is about 3 
times larger than for the attached dwelling units. 

Table 38 summarizes the mean home to work transport network 
distances for non-car owners and car owners for the fifteen Ontario 
census areas. This table shows that there are some large differences 
between the mean trip lengths for the two groups in some areas with the 
largest differences being in Hamilton, Sudbury, Ottawa, Sarnia and Windsor. 
The mean values presented in this table will differ from those presented 
earlier since all trips are included in the averages calculated in Table 38 
and not just the linkages between official census tracts. Some of the 
spatial differences in behaviour of the two groups have been presented 


in Chapter 5. 
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TABLE 38. Mean Trip Lengths for Non-Car Owners and 
Car Owners in Ontario Census Areas 


Census Area Mean Trip Length 
(km) 
Non-Car Car Owners 
Owners 
Guelph | 3.4 369 
Peterborough 330 3.4 
Sarnia 4.5 6.7 
Brantford 3.9 4.4 
Sault Ste. Marie 3%9 4.9 
Kingston 4.3 6.0 
Thunder Bay 4.6 5 .6 
Oshawa 369 4a 
Sudbury 6.4 ) 2 
Kitchener 4.1 bal 
Windsor 6.0 8.0 
London dal 629 
St. Catharines 4aJ 6.8 
Hamilton | 538 8.8 


Ottawa Gel 8.3 


= FAILS 


CHAPTER 7 


SOME GENERALIZATIONS ABOUT WORK TRIP MODELLING 


The overall objective of the studies of the 1971 census journey 
to work data described in reference [1] and in this report was to under- 
stand better the home to work spatial linkages patterns in the Ontario 
census areas. The basic motivation was the idea of developing general- 
ized relationships across 411 urban areas that could be used in transport 
planning studies. A traditional route has been followed in these studies 
consisting of the three components of (i) trip production analysis, (ii) 
trip attraction analysis, and (iii) trip distribution analysis. The 
hypothesis implicit in this approach was that the parameters produced 
by these three sub-analyses for each census area could be related to 
some broad properties of the census areas such as population size and 
employment sector composition. This chapter examines the extent to 
which these overall objectives have been met and outlines some of the 


implications for urban transport modelling. 


7.1 Trip Generation Analyses 


The trip generation analyses described in this report and the 
earlier study have demonstrated very clearly that high quality prediction 
equations for estimating the labour force composition of census tracts 
may be developed in terms of the dwelling unit composition of census 
tracts. The important feature of these prediction equations is that 


the partial regression coefficients are consistent across all of the 
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Ontario census areas with any real differences being due to variations 
in the labour force participation rates across census areas. These 
regression analyses have also shown that high quality prediction 
equations may also be developed in terms of census tract dwelling 
unit composition for estimating labour force transit captivity and 
non -captivity. Once again when differences in the labour force part- 
icipation rates between census areas are removed from the data the 
partial regression coefficients are very consistent across all census 
areas. 

These findings would suggest that work trip generation 
studies be standardized in terms of dwelling unit composition for 
all urban transport studies subsidized by the Ontario Ministry of 
Transportation and Communications. The dwelling unit composition 
and distribution in urban areas are perhaps the highest quality land 
use data available. More importantly it is a variable which is in- 
fluenced strongly by municipal planning policies and may be reliably 
estimated for future time horizons. Dwelling units are also the 
central item in turtieipal assessment records and the opportunities 
for developing more sophisticated dwelling unit classification schemes 
and relating travel behaviour to dwelling unit classes are enormous, 
Assessment records also contain information on household structure 
which would allow reliable correlations to be established between 
dwelling unit type and household structure variables such as labour 
force composition. 

Trip attraction equations cannot be established from the 


census data since information on the characteristics of places of 


— (oe 


employment is not collected. However, if a common system of employment 
and land use classification could be established across all Ontario 
municipalities then this trip attraction information could be exploited 
fully. An unsuccessful attempt was made to establish relationships 
between captive and non-captive trip attractions and employment sector 
type. Significant transport planning benefits would be derived from 

the adoption of a Province-wide municipal land use classification scheme 
and this may gradually emerge from the current studies of assessment 


practices, 


os 2 Trip Distribution Modelling 


The studies of trip distribution across the Ontario census 
areas have demonstrated that cross-sectional type gravity models are 
very poor in explaining the home to work patterns observed in the 1971 
census. A large number of modifications to the basic gravity model 
have been explored in this project and these have produced only marginal 
improvements to the explanatory power of the gravity model. However 
these analyses of the gravity model have provided important insights 
into its behaviour. 

The analyses which have been conducted have demonstrated very 
clearly that the best form of the gravity model to use is the doubly 
constrained version. However, the trip tables calculated by the gravity 
model when compared with the observed trip table exhibited goodness of 
fit levels equivalent to about a 75 percent randomly introduced error 
range. A detailed analysis of the calibrated doubly constrained models 


showed that the model properties are dominated almost completely by the 


Se 


spatial distributions of trip ends and that the so-called behavioural 
parameter of the deterrence function has only marginal effects. 

While the doubly constrained versions of the gravity model 
possess the best goodness of fit characteristics it is difficult to 
develop any rational interpretations of the production end and attraction 
end balancing factors or the deterrence function behavioural parameters 
either within or between census areas. Balancing factor magnitudes do 
not vary in any systematic way within urban areas and the deterrence 
parameter magnitudes have eee very weak association with mean trip 
length between urban areas. If the calibration parameters were rational 
then they should be interpretable in some way both within and between 
urban areas. 

Comparisons of the qualities of models calibrated in terms 
of network distances and network travel times showed that distance is 
the best measure to use. In very large urban areas the use of network 
times improved the goodness of fit of some of the outlying zones but 
the distance-based models were superior. Once again this observation 
reinforces the previous conclusion that doubly constrained model be- 
haviour is dominated by the structural parameters and not the behavioural 
parameter. 

The strategy adopted in this research project for improving 
the goodness of fit of gravity models was to ateeel ey trip ends and to 
distribute trip interchanges within each of these strata and then to 
sum the trip tables across strata. These attempts resulted only in 


marginal improvements to the doubly constrained models but there were 


- 205- 


some significant improvements to the production constrained models, 
Perhaps the greatest opportunity for improving the quality 
of gravity models is to stratify by period of development and to develop 
effectively a ent trip distribution model. The cluster analyses 
showed very clearly the importance of the Staging of development on 
the evolution of commuting sub-regions. This Strategy for improvement 
could not be implemented since information on the period of development 
of employment opportunities is not available from the census information. 
These studies of trip distribution across the fifteen Ontario 
census areas have demonstrated quite clearly that there is little point 
in calibrating and using the traditional cross-sectional type gravity 
models. The model-estimated base year trip tables are so different 
from the observed trip tables that any trip tables estimated for the 
horizon year are of little use. Time-staged or dynamic forms of the 
gravity model must be developed if estimates of future travel are to 
have any credibility. All data collection exercises should record the 
period of location at both the home and work ends in order that trip 
tables partitioned by period of development may be developed. It is 
also essential that every effort be made to collect journey to work 
data in the 1981 census in order to develop a consistent set of time 


series data on home to work linkages. 


- 206 - 


CONCLUSIONS AND RECOMMENDATIONS 


Consistent and high quality prediction equations for estimating 

the labour force composition of census tracts have been developed 
for all Ontario census areas in terms of the dwelling unit com- 
position of census tracts. 

Similar prediction equations have been developed for estimating 

the transit captivity and non-captivity of members of the residential 
labour force. 

Equivalent work trip attraction equations cannot be seston 
since independent measures of land use at the employment-end are 
not available. 

It is recommended that work trip generation studies be standardized 
in terms of dwelling unit composition. Centralized assessment 
records should be used to establish Ontario-wide dwelling unit 
classifications and consistent employment-end land use classifi- 
cations, 

Bi-proportional matrix balancing techniques may be used to establish 
the interaction structure with urban areas that is independent of 
zone size effects. Cluster analysis may be used very effectively 
to establish the commuting structure of urban areas and 1971 com-— 
muter-sheds have been established on a consistent basis for the 
Ontario census areas. 

Additional empirical analyses of other census data are required 

to establish more clearly the determinant of these commuter-sheds 


in the Ontario census areas. 
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The goodness of fit statistics such as the coefficient of deter- 
mination and chi-squared are inadequate for assessing the qualities 
of calibrated gravity models. The best statistic to use for com- 
paring observed and estimated trip matrices is the phi-statistic 
of information theory. 

Cross-sectional gravity models are very poor in explaining the 
observed journey to work patterns in Ontario census areas having 
goodness of fit characteristics similar to a trip table created by 
introducing random errors of 75 to 100 percent in the observed 
trip interchange magnitudes. 

Doubly constrained gravity models have Superior goodness of fit 
characteristics to the singly constrained models. With the doubly 
constrained models compliance with the trip end constraints is the 
major determinant of the estimated trip interchange magnitudes with 
the deterrence function having the secondary effect. The contri- 
bution of the deterrence function increases with urban area popu- 
lation. 

Gravity models stratified by geographic sub-region or by socio- 
economic group result only in marginal increases in the goodness 
of fit. The largest increase in goodness of fit was achieved by 
stratifying trip makers into car owners and non-car owners. 

The most important determinant of commuting that cannot be cap- 
tured through trip end stratification and constraint is the timing 
of development. Longitudinal models would be required to reflect 
timing of development effects in which trip interchange elements 


are partitioned into different sets to reflect different periods 
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of development. 

Network distances are adequate for use as a surrogate of general- 
ized travel costs in smaller and medium sized urban areas. Net- 
work travel times result in some improvements in the goodness of 
fit of gravity models in the largest cities. 

Because of the poor goodness of fit characteristics of the trad- 
itionally used gravity model it is recommended that their use as 
an estimator of future trip distribution patterms be discontinued 
in urban transport planning studies. It is recommended that the 
traditional cross-sectional model be replaced by a longitudinal 
form of the gravity model. Some development work is required in 


order to implement this recommendation. 
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